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ABSTRACT 


» 

This  research  is  a  study  in  seepage  analysis  for  engineered 
structures  built  below  the  groundwater  table;  in  particular  an 
analysis  of  underdrain  systems  for  such  structures.  General 
seepage  theories  and  design  considerations  are  discussed  as 
part  of  the  literature  review.  Also  design  concepts  and 
methods  of  analysis  are  covered  which  gives  a  brief  overview  of 
design  considerations  and  techniques  which  are  currently  used 
for  this  type  of  seepage  problem. 

Key  variables  to  the  design  of  an  underdrain  system  are  studied 
in  detail  by  using  a  two  dimensional,  steady  state  finite 
element  model.  Such  variables  studied  are  drain  spacing,  drain 
width,  anisotropic  hydraulic  conductivity,  thickness  of  the 
aquifer,  head  differential  and  sloping  aquifers. 

Results  are  presented  in  the  form  of  dimes: onl ess  plots  and 
nomographs  for  drain  discharge  and  free  surface  plots  for  the 
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CHAPTER  ONE 


INTRODUCTION 

Far  civil  engineered  structures  which  interact  with  sail,  a 
major  consideration  -for  design  is  the  presence  of  groundwater. 
Dams,  buildings,  dry  dock  -facilities,  underground  structures, 
pavements  and  earth  retaining  walls  must  be  designed 
considering  the  quantity  of  seepage,  the  position  of  the 
free  surface,  seepage  forces  and  seepage  control.  Good 
drainage  incorporated  in  the  design  of  these  and  other 
structures  can  often  benefit  through  substantial  cost 
reductions  as  well  as  increasing  the  safety  and  usefulness  of 
the  structure. 

Many  embarassing  and  sometimes  disastrous  problems  often  result 
from  man's  endeavors  to  control  seepage.  Saturated  soil 
backfills  have  caused  many  retaining  wall  failures  due  to 
increased  hydrostatic  pressures  following  extensive  rainfall . 
Highways  and  roadways  prematurely  deteriorate  and  fail  due  to 
trapped  water  in  the  pavement,  section.  Basement  slabs,  pond 
liners  and  drydock  facilities  have  been  severely  damaged  by 
hydrostatic  uplift  pressures  and  seepage.  Most  dam  and 
slope  failures  can  be  attributed  to  excessive  seepage  or 
uplift  pressures.  These  examples  and  untold  others  indicate 
that  good  seepage  and  drainage  control  together  with  an  adequat 
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understanding  of  groundwater  is  a  key  element  in  the  design  of 
most  civil  engineering  structures. 

Many  civil  engineering  structures  are  located  in  areas  where 
groundwater  causes  problems  of  uplift  and  seepage.  The  design 
therfore  must  consider  all  problems  associated  with  high 
groundwater  and  what  measures  must  be  taken  to  assure  the 
integrity  of  the  overal 1  design.  Casagr ande ( 1 965 ) ,  in  talking 
about  engineering  risks,  concludes  that,  "The  margin  of  safety 
that  we  incorporate  into  our  structures  should  bear  a  direct 
relationship  to  the  magnitude  of  potential  losses,  and  it  must 
also  take  into  account  the  range  of  uncertainty  involved." 

This  important  concept  hopefully  will  be  considered  throughout 
the  presentation  of  facts  and  data  obtained  in  this  paper. 
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Prior  to  the  20th  century,  analysis  of  seepage  and  the  control 
of  drainage  for  enginerring  structures  was  done  by  experienced 
engineers  who  based  their  designs  on  trial  and  error  methods 
which  sometimes  were  successful  and  sometimes  not.  Though 
experience  has  proven  to  be  very  neccesary  and  helpful  in  the 
design  of  drainage  systems,  very  useful  tools  employed  today 
are  numerical  methods  such  as  finite  element  and  finite 
difference  techniques  which  require  the  use  of  the  digital 
computer.  Their  problem  solving  capabilities  have  made 
possible  solving  very  complicated  and  complex  groundwater 
pr  ob 1 ems . 
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The  purpose  of  this  research  is  to  study  the  design  of 
underdrain  systems  for  large  structures  placed  below  the 
groundwater  table.  In  chapter  two,  basic  seepage  theories  and 
design  requirements,  designs  concepts  and  useful  methods  of 
analysis  are  presented  as  a  brief  literature  review.  This 
chapter  covers  critical  aspects  of  designs  that  should  be 
considered  and  contains  a  brief  review  of  methods  used  by 
geotechnical  engineers. 

In  chapter  three,  the  finite  element  program  developed  by  Dr. 
M.  Aral,  of  the  School  of  Civil  Engineering  at  Georgia 
Institute  of  Technology,  is  discussed  describing  how  the 
design  of  building  underdrain  systems  can  be  modeled  using 
this  technique.  This  finite  element  program  made  possible 
studying  effects  of  key  variables  such  as  permeabi 1 i ty, 
hydraulic  head,  drain  geometry,  depth  of  groundwater,  and 
anisotropic  soils.  Program  verification  is  performed  using 
another  program,  SEEP84,  modified  by  the  Civil  Engineering 
Department  at  Virginia  Polytechnic  Institute  and  State 
University,  for  the  microcomputer.  Dupuit’s  solution  is  also 
used  to  compare  selected  cases  with  the  finite  element 
progr ams . 


Chapter  -four  presents  data  generated  -from  the  -finite  element 
study  such  as  position  o-f  the  -free  surface  and  -flow  quantity 
■for  selected  systems  of  underdrains  suitable  for  use  beneath  a 
foundation-  The  data  is  presented  in  the  form  of  charts  and 
nomographs  to  be  used  as  a  guide  for  underdrain  system  design. 

Chapter  five  summarizes  the  results  of  the  study  and  presents 
appropriate  concl ussi ons. 

Appendix  A  contains  program  documentation  for  the  finite 
element  program  used  in  this  study  which  details  to  the  user, 
documentation  to  the  use  of  the  program.  Appendix  B  presents 
example  input  and  output  data  generated  by  this  study. 
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CHAPTER  2 


LITERATURE  REVIEW 


2.1  Seepage  Theories  and  Design  Requirements 


Groundwater 


Soils,  sometimes  described  as  a  porous  media,  is  a  mass  made  up 
of  solid  matter  having  interconnected  voids,  passages,  -fissures 
or  pore  spaces.  These  pore  spaces  can  be  filled  by  a  liquid 
and  or  a  gaseous  phase.  Figure  2.1  shows  several  typical  soil 
and  rock  indices  which  depicts  how  pore  spaces  can  act  as 
conduits  for  water  or  subcapillary  openings  in  which  water  may 
be  held  by  adhesive  capillary  forces  < Mei nz er , 1 942 ) . 

The  soil  mass  forms  a  complex  network  of  irregular  passages 
through  which  the  fluid  may  flow.  A  key  parameter  of 
groundwater  movement  is  the  porosity  of  the  soil.  As  gravity 
acts  on  the  fluid,  it  moves  through  the  pore  spaces  at  varving 
rates  through  each  open  pore  space.  From  a  macro  scale,  this 
flow  has  been  mathematically  described  by  several  authors 
(Bear,  1981;  Terzaghi,  1943;  Cedergren,  1977;  Harr,  1962; 

Bouwer ,  1978; Freeze,  1971). 


Figure  2.1  Diagram  showing  soil  and  rock  interstices.  A.  We 
sorted  sand  deposit  having  high  porosity;  B.  Poorly  sorted 
deposit  having  low  porosity;  C.  Particles  that  may  be  porous 
themselves;  D.  Deposit  whose  porosity  may  be  diminished  by 
deposition  o-f  mineral  matter  in  the  interstices;  E.  Rock  that 
is  porous  by  solution;  F.  Rock  that  is  porous  by  fracturing 
(after  Meincer,  1942) 


Moisture  content  in  soils  can  vary  with  depth  -from  the  surf ac 
as  shown  by  Figure  2. 2.  Most  seepage  problems  considers  the 
zone  o-f  saturation  which  is  the  zone  where  the  soil  mass  is 
completely  saturated.  The  terms  groundwater  level,  phreatic 
water  level,  water  table,  -free  surface  or  groundwater  surface 
are  all  synonymous  to  denote  the  level  o-f  water  where 
atmospheric  pressure  exists. 
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In  dealing  with  capillary  -fringe  in  fine  grained  soils  the 
height  of  water  can  very  significant.  As  the  pore  diameter 
decreases,  the  capillary  rise  increases.  Although  the  soil  l 
considered  unsaturated,  the  soils  moisture  content  will  be 
higher  with  depth. 

The  cone  of  capillary  rise  and  fringe  zone  may  contain  large 
quantities  of  water  for  fine  grained  soils.  Water  removal  by- 
drainage,  due  to  capillary  forces,  may  be  difficult  for  such 
fine  grained  soils.  Figure  2.3  shows  how  the  actual  moisture 
distribution  is  appro::  i  mated  by  a  step  distribution.  This 
appr ox i mat i on  of  the  capillary  fringe  level  above  the  free 
surface,  hc  ,  is  often  neglected  in  some  seepage  problems. 
Several  emperical  methods,  based  on  grain  size  diameter  and 
porosity  of  soil,  are  available  to  estimate  hc  (Mavis  and 
Tsui,  1939;  Pol ubari nova,  1952,  1962;  Silin  Bekchurin,  1958). 
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Figure  2.3  Appro::  i  mat  i  on  of  phreatic  surface  and  capillary 
-fringe  zone.  (after  Bear,  1979) 

Many  textbooks  on  seepage  and  groundwater  discuss  water  bearing 
and  retention  capibilities  o-f  the  capillary  -fringe  and 
movement  of  moisture  through  this  zone  (Terzaghi,  1943;  Kovacs, 
1981;  Sowers,  1979).  However,  when  this  height  of  water  is  very 
small  relative  to  the  total  thickness  of  the  aquifer,  the 
capillary  fringe  is  sometimes  neglected.  The  geotechnical 
engineer  should  be  concerned  with  the  unsatur ated-satur ated 
flow  condition  when  the  stability  of  the  engineered  structure 
is  to  be  analyzed,  such  as  a  dam. 

Groundwater  analysis  for  seepage  problems  must  consider  the 
quantity  of  water  to  be  drained,  the  soils  resistance  to 
drainage,  and  what  effects  the  drainage  may  have  on  the  soil 
and  the  structure.  Movement  of  groundwater  can  be  categor  l  ze-:1 


9 


either  as  seepage;  capillary  rise  or  flow;  percolation;  or 
turbulent  flow.  For  this  study,  only  movement  of  groundwater 
by  seepage  is  considered. 

Design  Requirements 

Structures  built  below  the  watertable  must  be  designed 
considering  the  following: 

(1)  structural  ability  to  withstand  the  applied  hydrostatic 
water  pressures, 

(2)  drainage  control  which  short-cuts  the  flow  or  movement 
of  groundwater  from  getting  to  the  structure, 

(3)  drainage  control  which  is  incorporated  with  the  overall 
design  to  remove  the  groundwater, 

(4)  assure  water  velocities  be  small  enough  to  prevent 
excess  particle  migration  and  erosion. 

Savings  of  dollars,  time  and  natural  resources  can  be 
recognized  when  structures  can  be  considered  on  dry,  stable 
material ; in  addition  to  thinner  basement  slab  thicknesses  may 


In  1856,  Henry  Darcy  investigated  the  -flow  of  water  in  vertical 


homogeneous  sand.  This  work  has  developed  the  most  -fundamental 
relationship  which  is  still  used  in  seepage  analysis  today. 
Darcy’s  Law  states: 


Q  =  k  i  A 


(2.  1 ) 


where  Q  (L  Vd  is  a  volume  discharge  or  flow  of  a  given 

2 

cross-sectional  area,  A  <L  )  in  a  given  time.  The  discharge  is 
also  a  function  of  the  hydraulic  gradient,  i  (L/L)  and  a 
resistance  coefficient,  k  <L/T),  the  permeability  of  the  soil. 
Investi gators,  Muskat  (1937),  Taylor  (1949)  and  Leonards 
(1962),  have  concluded  that  Darcy’s  Law  is  valid  where 
velocities  are  low  or  when  flow  is  laminar,  i.e.  Reynolds 
number  less  than  10. 


Hydraulic  Conductivity 


The  soils  resistance  to  flow  or  the  coefficient  of 
porporti onal i ty ,  k,  for  Darcy’s  Law,  is  also  called  the 
hydraulic  conductivity.  Hydraulic  conductivity  is  a 
coefficient  which  expresses  the  ease  or  resistance  the  water 
has  flowing  through  the  soil.  It  depends  on  the  soils  matrix 
and  fluid  properties  of  the  permeant. 


Many  researchers  have  performed  tests  to  attempt  to  develops 
empirical  rel at i onsh 1 ps  of  the  soils  properties  to  hydraulic 
conductivity.  Some  of  these  properties  considered  are  grain  or 
pore  size  distribution,  shape  of  the  grains  or  pores, 
tortuosity,  specific  surface  and  porosity  (Mitchell,  1976s  K'ozeny 
1927;  Carman,  1956;  Cedergren,  1977;  Kovacs,  1981).  All  report 
that  careful  field  i nvesti agati on  and  careful  conductivity 
studies  are  necessary  for  each  design.  Figure  2.4  is  presented 
to  show  some  characteristics  of  various  soil  types  and  most 
important,  methods  for  determination  of  selected  properties 
(Terzaghi  and  Peck,  1948). 

In  the  unsaturated  zone,  the  hydraulic  conductivity,  moisture 
content  and  pressure  head  is  extremely  complex  and  not  well 
understood.  Again,  field  measurement  is  the  best  method  for 
obtaining  the  best  parameters,  such  as  the  use  of  theistor 
psychrometer  (Kay  and  Low,  1970)  or  gamma  ray  absorption 
techniques  (Ferguson,  1970;  Aral ,  1983).  However,  many 

empirical  formulas  have  been  developed  to  relate  the  soil 
parameters  to  a  relative  hydraulic  conductivity  of  the 
un-saturated  zone  (Kovacs,  1981;  Irmay,  1954;  Childs  and 
Col  1  i  s-Geor  ge ,  1950;  Gardner,  1958;  Brooks  and  Corey,  i9m>; 

Maslia,  1980;  Reeves  and  Duriguid,  1975). 


2.2  Design  Concepts 


Seepage  Forces 

The  effects  of  water  on  the  soil  and  the  structure  is 
controlled  by  the  pressures  in  the  water.  In  soil  mechanics, 
the  water  pressure  is  defined  as  neutral  stress.  For  water  not 
moving,  or  steady  state,  the  pressure  can  be  computed  from  the 
basic  equation  of  hydrostatics 

h  =  P  +2  (2.2) 


where  h,  energy  head,  is  equal  to  the  positic  ,  2,  plus  the 
pressure  head  or  piezometric  head,  P/^  .  Typically  in 

seepage  designs,  the  velocity  head  is  very  small  and  mav  be 
neglected.  Then  the  pressure  head  plus  elevation  is  equal  to 

head,  h  (Cedergren,  1977). 

For  the  condition  that  water  pressure  is  greater  than  the  total 
load,  say  of  the  soil  and  structure,  the  pressures  are  termed 
uplift  pressures  and  liquifaction  can  occur.  Therfore  the 
control  of  uplift  pressures  must  be  a  primary  design 


con si der at i on . 


Critical  Velocity 

According  to  Darcy's  Law,  the  discharge  velocity  o-f  seepage 
through  soi 1  is 

v  =  k  x  i  -  (2.3) 

d 

Seepage  is  only  occuri ng  through  the  pores  and  voids. 

Therefore  the  true  seepage  velocity  should  be  a  -function  of  th 
porosity  of  the  soil  mass.  Re-writing  equation  (2.3)  the 
seepage  velocity  can  be  expressed  as 

k  x  i 

v  =  -  (2.4) 

sn  ^ 

e 

where  r\  e  is  the  effective  porosity. 

In  flow  conditions,  the  saturation  of  soil  is  moving  into  or 
away  from  the  unsaturated  soil.  The  capillary  head.  hc  ,  is 
worfing  with  gravity  to  increase  the  hydraulic  gradient.  For 
this  condition  the  hydraulic  gradient  can  be  expressed 

h  +  h 

i  =  _ 9 

L  (2.5) 


Fig  u r e  2.5 


depicts  this  condition  -for  -flow  through  a  sand. 


Reservoir 


Reservoir 


Screen 


Figure  2.5  Force  of  saturation. 


Boiling  occurs  when  seepage  pressures  in  an  upward  direction 

exceeds  the  downward  -force  o-f  the  soil.  This  must  be  analysed 

•for  all  seepage  designs.  Critical  hydraulic  gradient,  i  ,  can 

c 

be  expressed 


where 


*w 

*b 


i  c 


(2.6) 


total  unit  weight  of  soil 
unit  weight  of  water 

buoyant  unit  weight  of  soi 1 ( NAVFAC , DM -7 . 1 ) . 
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Fi 1  ter s 


A  very  important  design  consideration  is  the  possibility  of 
seepage  -forces  moving  erodible  soils  and  rock,  termed  piping. 
Betram  (1940),  under  Terzaghi  and  Casagrande,  established  the 
■following  criteria  for  -filter  design,  based  on  grain  size. 


D 


D 


15 


85 


( of 


(of 


filter) 


so  i  1  ) 


<  4 


to  5 


< 


D  ^  (of  filter) 

D  (of  soi 1 )  (2.7) 


This  relationship  is  represented  by  Figure  2.6  to  illustrate 
the  prevention  of  piping  by  filters  (Cedergren,  1977). 


LEGEND 


I  =  in-piace  soil 


U,-,  (so  l) 


=  D-\  soil  particle, 
entrapped  in  filler 

=  soil  which  has 
-migrated  into  filter 
and  is  held  by  £),.-, 
size  soil  particles 


\ 


\ 


\ 


\  y 


\ 


I  -  . 

*  r:' 


Soil  > 


Nom.nal  boundary 
before  stabilization  - 
ur.cier  seepage 


‘  o 


/,■ 


■ — 
.  / 


,  L)\-  C.iter) 

) 


Filter 


Figure  2.6  Illustration  of  prevention  o+  piping  by  filter-.;, 
(after  Rankilor,  1981) 
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Geosynthet ics 


Today,  new  synthetic  membranes  are  used  to  serve  the  purpose  of 
either  filtration  or  separation.  This  two  factors  are  very 
important  parameters  for  most  seepage  problems.  Figure  2.8 
shows  comparisons  of  some  of  the  features  of  membrane-wrapped 
drains  to  conventional  granular  pipe  drains.  Also  Figure  2.9 
shows  the  bridging  effect  caused  by  good  seperation  of  the 
larger  materials  and  smaller  particles. 


firannlar  drain 


Membrane-wrapped  drain 


*  +  9  ^ 


Pipe  drains 
need  frequent 
rod  ding  points. 
Expensive. 


Pipe  must  be 
levelled  by 
engineer,  to 
be  an  exact 

gradient.  This  is  expen¬ 
sive.  Trench  must  be  left 
open  for  this,  and  can  get 
silted  up  during  storms. 

Well-point 
pumping 
may  be 
necessary  to 
lay  the  pipe. 


Sand  and  stone  must  be 
supplied  and  laid  to  a  strict 
specification.  Very  expensive. 

Pipe  Is  expensive  to  buy,  and 
must  be  laid  by  hand.  Trench 
must  therefore  be  wide,  and 
timbered  if  deep. 


a> 

o 

>_  o. 
a)  » 

II 


Wide  trench 
means  more 
excavation  and 
more  backfill  than  Is 
technically  necessary 


The  deep  transporter  means 
that  settlement  can  take  place  ' 
without  the  drain  blocking  — 
unlike  a  pipe 

Can  be 


taken  deeper 
because  sides 
do  not  need 
support 


▼ 


±_r_ 


Cos'  r>-,  shoring 
up  open  trench 
is  eliminated 

Granular  fill  only 
supports  trench, 
so  can  be  cheap 
specification 

This  drain  is  very 
much  quicker  to 
construct  &  can 
be  laid  in  flooded 
ground. 

Narrower 

width  because  no  pipe  to  be 
laid,  so  men  do  not  have  to 
go  into  the  trench. 


Figure  2.8  Conventional  granular  pipe  trench  drain  and 
membr ane-wr apped  trench  drain.  (after  Rankilor,  1981) 
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A  goad  representation  of  how  geosynthetics  combine  action  of 
seperation  and  filtration,  is  shown  by  Figure  2.10.  As  water 
has  more  area  of  flow  with  the  membrane,  the  seepage  can  be 
considered  to  be  at  a  faster  rate  or  drain  quicker  (Rankilor, 
1981 ) 


Without  membrane 


With  membrane 


i 


/  Pressure  .  \  ♦  L L^v/ 

<  /  T  v T  '  Pressure  T  rf 

I  I  1  '(  . 

Water  tries  to  \ — ' 

escape  through  «.crTinid  Clay 

small  area  SECTION  ,oj) 


Water 


escapes  SECTION 


rapidly 


o  o  “  o  o 

Clay 

O  o  o  so,'c 

k  0/ 


Without  membrane 


Small  surface  area  Large  stone/clay 

for  moisture  loss  contact  area 


With  membrane 


Large  surface  area 
for  moisture  loss 


Small 
stone/clay 
contact  area 


Figure  2.10  Seperation  effects  by  using  a  membrane, 
(after  Rankilor,  1981) 


i  ^  "  «  ’  «  *  ■  *  «_*  «*,. 


2.3  Methods  of  Analysis 

In  analyzing  seepage  problems;  geometry,  boundary  conditions, 
soil  parameters,  flow  quantities  and  location  of  the  free 
surface  must  be  quantified  by  some  means.  Though  this  involves 
many  considerations,  the  engineer  must  gain  reasonable  safe 
results  for  the  design. 

Through  the  history  of  seepage  analysis  many  techniques  and 
methods  have  been  developed,  such  as  graphical,  analytical  and 
numerical  solutions.  Still  used  today  is  the  graphical  method 
of  flow  net  construct i on .  However  with  the  complex  geometries 
envolved  and  via  the  speed  of  computers,  numerical  models  have 
been  developed  as  very  useful  tools  for  the  engineer. 

Flow  Nets 

Many  researchers  have  developed  methods  for  analyzing  flow 
through  parous  media.  The  oldest,  developed  by  Caaagrande  and 
Terzaghi  (1937),  is  the  use  of  flow  nets,  a  graphical  approach 
to  seepage  analysis. 
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The  continuity  equation  in  the  -form 


3u  3v  3w  (2.7) 

-  +  -  +  —  ~  0 

3x  3y  3c 


shows  {or  a  Cartesian  coordinate  system,  the  quantity  of  water 
entering  an  element  must  equal  that  leaving.  The  terms  u,  v 
and  w  are  discharge  velocity  components  o-f  directions  a,  y  and 
z •  respectfully. 


According  to 

Darcy’s  Law,  expressed 

by  equation  <  2 . " 

3)  ,  the 

component  of 

discharge  velocity 

can 

be  expressed 

3h 

3h 

3h 

u  =  - 

-k  -  ,  v  =  -k 

— 

,  w  =  —  k  - 

(2.8) 

3  m 

3y 
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Subst i tut i ng 

into  equation  (2.7) 

—  k  <  3  h  /  3 : 

•;)  -k(3h/3y) 

-k 

( 3  h  /  3  z  ) 

(2.9) 

a  - 

-  +  a - -  + 

a  - 

-  =  0 

3m  3  y  3  z 
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if  m  .1  k  I\  • 


If  k  is  constant  -For  the  element  and  then  expressed  in  a 


typical  two-dimensional  form,  equation  (2.9)  becomes  the  common 
Laplace  equation. 


a2h 


92h 


(2.10) 


~2  + 


=  0 


9y 


The  general  conditon  of  flow  analysis  using  flow  nets  is 
represented  by  Figure  2.11.  Rules  and  good  examples  are 


presented  in  many  textbooks  (Cedergren,  1977;  Harr,  1967). 


free  water  surface 


Figure  2.11  General  condition  for  flow  lines  and 
equi potent l al  lines  in  flow  net  construct i on . 
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Dupm  t  Assumption1 
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For  unconhned  seepage  problems  the  engineer  must  calculate 
the  location  of  the  -free  surface.  Many  problems  of  flow  have 
been  appr a:  1  mated  based  on  the  Dupuit  theory  of  unconf ined 
flow.  Dupuit  (IBF’3 )  assumed;  '1)  for  small  inclinations  of 
the  line  of  seepage,  the  streamlines  can  be  t.  afen  as  horizontal 
and  equal  potential  1  ines  vertical,  arid  (2)  the  hydraulic 
gradient  was  equal  to  the  slope  of  the  free  surface  and 
invarient  with  depth  (Harr-,  1962). 


In  the  case  of  Figure  2.12  the  Dupuit  -formula  is  expressed  as 


2L 


Many  other  reports  and  textbooks  have  presented  charts  and 
nomagraphs  based  on  these  assumptions.  Later  in  this  study  a 
nomograph  -for  a  simple  trench  drain  with  blanket  drain 
solution  developed  by  the  Federal  Highway  Admi n 1 strat i on  <1980) 
will  be  presented. 

Numerical  Analysis 

Flow  through  parous  media  has  been  appr ox i mated  by  many  authors 
using  ordinary  non-linear  differential  equations  and  partial 
differential  equations.  Coupled  with  the  speed  of  the  computor 
numerical  solutions  to  many  seepage  problems  have  been 
developed.  It  should  be  pointed  out  that  complex  mathematical 
problems  solved  by  numerical  methods,  are  only  as  accurate  as 
the  math  model  can  represent  the  conceptions  and  phenomina  of 
seepage  flow.  In  most  cases,  simplification  is  the  rule,  and 
numerical  models  can  only  be  treated  as  approximations  to  the 
actual  occurances.  Nevertheless,  the  number  crunching 
capibilities  and  time  savings  involved  gives  the  engineer  much 
needed  time  for  the  true  analysis  of  the  seepage  problem. 

The  analytical  and  graphical  methods  mentioned  previously  are 
sometimes  difficult;  and  limited  to  flow  systems  which  are 


relativly  simple.  In  actual  -field  conditions,  the  soils  may  be 
heterogeneous  and  anisotropic  and  boundary  conditions  very 
complex.  These  conditions  make  closed  -form  methods  very 
difficult  and  many  simplifications  would  have  to  be  made. 
However,  with  the  numerical  models,  these  complexities  can  be 
easily  handled  and  many  more  check  cases  can  be  made  in  much 
less  time. 

Finite  Difference  Methods 

Many  numerical  computer  models  have  been  developed  using  finite 
difference  methodology.  Since  the  classic  work  of  F'eaceman  and 
Rachford  (1955),  the  oil  industry  has  sponsored  many  studies  on 
model  development.  A  good  historical  background  of  the  finite 
difference  methods  developed  in  seepage  analysis  is  giving  by 
Huyakorn  (1983).  Worth  mentioning  are  the  works  of  Cooley 
(1971)  and  Freeze  (1971)  in  the  development  of  other  finite 
difference  models.  These  models  solve  for  the  free  surface  and 
treat  the  saturated-unsatur ated  flow  as  an  immobile  air  phase. 

Finite  Element 

The  finite  element  method  for  salving  subsurface  flow  has  been 
widely  used  and  researched  be  soils  engineers  and  groundwater 
hydrologists.  The  program  used  in  this  research  and  discussed 
in  more  detail  in  Chapter  3  is  a  finite  element  model. 


host  -finite  element  models  are  built  on  the  Galerkin  method. 


* 
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Some  earlier  -finite  element  models  -for  steady  state  and 
saturated  flow  were  developed  by  Desai  and  Abel  (1972),  Conner 
and  Brebbia  (1977)  and  Zienkiewicz  (1977).  Finder  and  Frind 
(1972)  used  isoparametric  quadr i 1 ater al  elements  with  the 
llalerkin's  technique.  Neuman  (1973,  1974)  used  a  tour  node 

quadr l 1 ater al  element.  A  good  reference  for  other  works  in  the 
finite  element  models  that  have  been  developed  is  reported  in 
Huyakorn  and  Finder  (1983)  and  Heijde  (1985). 

Some  significant  works  were  tho*«  of  Taylor  and  Brown  (1967) 
where  a  finite  element  model  was  developed  for  seepage  through 
a  dam  solving  for  the  free  surface.  The  flow  was  only 
considerd  for  the  saturated  zone  *n  this  model.  Freeze  (1971) 
demonstrated  that  the  conditions  of  the  unsaturated  zone  may 
strongly  influence  the  position  of  the  free  surface  and  results 
obtained  not  considering  the  capillary  action  for  this  case 
could  be  in  error.  Fredlund  and  Morgenstern  (  1976)  devel oped 
improved  constitutive  equations  for  an  unsaturated  soil.  Lam 
and  Fredlund,  etal  (1984)  presented  additional  theories  and 
numerical  development  for  further  modeling  techniques  to 
consider  the  unsaturated  zone  in  seepage  analysis. 


CHAPTER 


FINITE  ELEMENT  STUDY 


;.l  Finite  Element  Program  Theory 


Introducti on 


Finite  element  analysis  of  seepage  problems  can  be  described  as 
follows  (Huyakon  and  Finder,  19B3) : 

(1)  A  discrete  number  of  i nter connected  nodes  to  define  a 
series  of  sub-regions  or  elements.  These  can  take  the  form  of 
triangles,  quadr i 1 ater al  or  rectangular  shapes  which  are  so 
spaced  to  define  the  geometric  boundaries  and  soil  types. 

(2)  A  matrix  expression  is  developed  to  relate  flow 
variables  to  each  node  of  the  elements.  The  mathematical 
expression  may  be  of  the  form  of  either  a  variational  or 
weighted  function.  A  weighted  residual  function  is  used  in  the 
Galerkin  method  (Conner  and  Brebbia,  1977). 

(3)  A  set  of  algebraic  equations  that  describe  the  entire 
global  system  is  applied  to  assemble  the  element  matrices. 

These  equations  are  derived  to  consider  compatible  conditions 
at  each  node  shared  be  each  adjoining  element. 

(4)  Boundary  conditions  of  flow  are  incorporated  to  reflect 
appropriate  limits  to  the  global  matrix  equations. 

(5)  The  resulting  set  of  simultaneous  algebraic  equations 
are  solved  by  some  alogrithm  method.  Many  algor i thyms  have 


i/.  Pi 


» 
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been  used  to  efficiently  take  advantage  of  the  banded  and 
symmetric  -features  of  the  coe-f -f  i  c  i  ent  matrix  (Bathe  and  Wilson, 
1976)  . 


The  -finite  element  method  has  great  potential  in  seepage 
analysis  problems,  in  that  the  complex  geometries  o-f  flow 
regions  can  be  easily  constructed  via  nodes  and  elements.  Th 1  ■ 
chapter  will  present  the  general  theories  and  formulation  of 
equations  for  the  model  used  in  this  study. 


Program  Theory 


A  finite  element  program  developed  by  Dr.  Mustusfa  Aral  of  the 
School  of  Civil  Engineering,  Georgia  Institute  of  Technology, 
was  used  to  study  the  flow  and  variables  involved  for 
subsurface  drains  as  earlier  described.  This  program  was  used 
to  see  effects  of  different  parameters  in  the  analysis  performed. 
The  program  was  operated  on  the  Georgia  Institute  of  Technology 
main-frame  CDC  CYBER  computer.  A  complete  user’s  guide  for  the 
program  is  included  as  Appendix  A-l.  Also  example  input  and  output 
used  in  this  study  is  included  by  Appendix  B-l.  The  formulation 
for  the  program  presented  here  was  summarized  from  a  previous  stud'. 
(Aral  ,  Sturm  and  Fulford,  1981). 


This  model  only  considers  the  flow  condition  of  the  free 
surface  being  the  point  where  the  water  level  is  at  atmospheric 
pressure.  In  other  words,  this  model  does  not.  treat  the 


unsaturated  zone  o-f  the  soil. 


For  this  study  the  -flow  continum 


M 


! 

$ 
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is  modeled  to  determine  the  quantity  o-f  -flow  at  each  drain  and 
to  determine  the  -free  surface  location  between  each  o-f  the 
drains.  This  model  also  computed  the  velocities  at  each 
element  -for  a  good  check. 

Equations  of  Flow 

Darcy’s  Law  is  used  as  the  basic  relationship  governing 
groundwater  -flow.  Flow  is  assumed  to  be  laminar  and  the 
inertial  -farces  are  negligible  compared  to  viscous  -forces. 
Rewriting  Darcy’s  Law  in  the  cartesian  coordinate  system: 

-*■  -*■ 

v  =  -  K  V  h  (3.1) 


1 


V* 


where  V  is  seepage  velocity,  h  is  piezometric  head  and  K  the 
hydraulic  conductivity.  In  the  two  dimensional  case,  equation 
(3.1)  can  be  written  as 


9  h 

v  -  -  Ki±  -  ,  i-  1,2  (3.2) 

9  x 


where  is  a  conductivity  tensor. 
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The  continuity  equation  -for  the  -flow  in  an  element  of  soil  can 
be  expressed 

9  9 

_  <°V.)  =  - _ <n  D  >  ,  1  =  1,2  (3.3) 

3>;  3t 


where  h  is  porosity  of  the  soil  and  p  is  the  density  of  the 
•fluid.  Assuming  the  -fluid  and  soil  particles  are 
incompressible,  equation  (3.2)  and  (3.3)  can  be  used  for  the 
following  expression  in  the  two-dimensional  form: 


9  3  h 


Figure  3.1  can  be  refered  to  for  notation  of  the  key  variables 


used  is  this  study. 


Boundary  Conditions 


For  -free  surface  problems,  which  is  the  case  for  this  study, 
the  enact  location  of  the  phreatic  surface  is  not  known  at  the 
beginning  of  the  solution.  Location  of  the  free  surface  is  an 
iterative  process  using  two  independent  boundary  conditions 


available  at  the  free  surface.  The  line  of  seepage  at  the  free 
surface  is  the  upper  streamline  of  the  flow  domain.  At  every 


point  along  this  streamline  the  pore  pressure  is  at  atmospheric 
pressure  and  theref ore  the  piezometric  head  equals  the  elevation 
head.  Also  the  flow  across  this  boundary  is  zero  because  of  it 
being  a  streamline. 


Considering  the  seepage  domain  as  depicted  by  Figure  3.1,  the 
steady  state  boundary  conditions  can  be  described  as  follows: 

(a)  Impervious  Boundary  (AB) :  The  program  considers  this 

boundary  as  a  streamline  where  velocity  to  this  boundary  equals 


(b)  Constant  Head  Boundaries  (AC)  and  (EF):  Constant 
hydrostatic  pressure  is  exerted  on  these  boundaries  where  total 
piezometric  head  is  constant,  i.e.  equal  to  the  water  level  in 
a  piezometric  tube. 

<c)  Phreatic  Surface  (CE) ,  (FG)  ,  (HI),  (JK)  and  (LM) :  This 

boundary  is  a  combination  of  both  cases  (a)  and  (b)  above.  By 
iteration,  described  later,  both  conditions  must  be  satisfied 
independent! y. 


( d )  Seepage  Faces  (EF),  (GH),  (IJ),  (KL)  and  (MN> :  Seepage 

forces  occur  at  these  boundary  faces  where  water  exits  through 
the  soil  continue,  i.e.  trench  drains.  Therfore  streamlines  of 
flow  must  be  normal  to  these  boundaries. 
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3.2  Finite  Element  Formulation 

General  Equations 

The  -first  step  in  solving  the  seepage  problem  by  the  -finite 
element  method  is  to  discretize  the  -flow  region  into  a  finite 
number  of  subregions  called  elements.  Using  finite  elements, 
equilibrium  equations  must  be  made  to  handle  the  geometries 
involved.  This  results  in  a  set  of  simultaneous  algebraic 
equations  rather  than  differential  equations.  Figure  3.2 
depicts  the  quadri lateral  elements  used  for  this  program  and 
shows  how  four  six  nodal  triangles  are  formed  internal  to  the 
program.  A  more  detailed  description  is  of  finite  element 
equations  and  formulations  can  be  found  in  Zienkiewicz  (1971) 
Desai  (1972)  and  Aral  (1974,  1976). 


Figure 


Typical  element  and  nodes 


A  finite  element  approximation  to  equation  <o .  5 )  can  be  written 
through  a  variational  approach.  Assuming  Dirichlet  and  Neuman 
boundary  conditions,  and  considering  the  two  dimensional  case 
of  steady  flow,  the  fallowing  equations  can  be  expressed: 


on  S  , 


(3.  7) 


and  Vn  -  IK  u  ^  anl  +  K22  ^  “nd  _  Vn  on  S2  (3-7) 


where  h  is  the  specified  Direchlet  boundary  condition  in  terms 
of  piezometric  head;  VR  is  the  velocity  in  the  direction 
normal  to  the  boundary;  and  cx^  are  the  direction  cosines  of 

the  normal  with  respect  to  and  ^  ?  and  and  S2  indicate 

respective  boundaries  of  the  domain.  From  equations  (3.5)  and 
(3.7),  the  fallowing  variational  statement  can  be  written: 


{hd„id>!2 


A  f  3-, . 


[v-  -  V»1 


<5  hdS 


(3.8) 


Integrating  by  parts  one  obtains 


{  9h  3  6h  3  h  3  6h 

)KII  +  2  2  ' 

A  (  *1  9*1  9  x2  9  Hi 


V  6  hdS  =  0 

n 


dx  i  d;;  2 


(3.  9) 


Since  -from  earlier  discussed  boundary  conditions  as  -for  this 
problem  studied,  Vn  =  0.  Thus,  equation  (3.9)  can  be  further 
simplified  by  eliminating  the  boundary  intregrals. 


Assuming  equation  (3.9)  is  written  for  a  single  element,  matri; 
equations  can  be  generated  using  an  i nterpol at l on  function  for 
appro::  imati  ng  h  over  an  element.  The  interpolation  function 
used  in  this  model  is  given  as 


h  =  Nk  hk 


k  =  1 , 2,3, . . .6  <3.10) 


where  h^’s  are  the  unknown  nodal  values  of  the  dependent 
variable  and  N ^  is  the  i nterpol at l ng  polynomial  used  to 
approximate  h  in  an  element. 


For  a  typical  si  :-nodal  triangle  element  these  polynomials  are 


g 1 v en  as 

N  t  -  ?i<2Ci  -1)  i  =  1,2,3  (3.11) 

and  N  .  +  3  =  4^  j  =2,3,1  (3.12) 

where  ^  ,  and  C3  are  area  cordinates  defined  by  the 

following  relations 


where  A^  ,  A2  and  A  3  are  the  areas  of  the  three  subtr 1 angl es 
subtended  by  the  point  (F')  and  the  corners,  the  index  on  (A) 
designating  the  opposite  corner  numbers,  and  A  is  the  area  of 
the  triangle  (123)  shown  in  Figure  3.3. 
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To  simplify  the  resultant  equations,  the  origin  of  the  local 
axis  which  defines  the  local  coordinates  of  the  element,  nodes 
are  placed  at  the  centroid  of  the  element  and  the  principle 
axis  are  inclined  angle  @  in  the  direction  of  local 
an i sotropy . 


Substituting  equation  (3.10)  into  equation  (3.9)  the  final  form 
of  the  approximating  equations  in  an  element  can  be  expressed. 
The  stationary  condition  of  equation  (3.9)  is  obtained  by 
equating  the  first  partial  derivatives  of  the  function,  with 
respect  to  the  undetermined  parameters,  hj  ,  to  zero: 


3N*  3Nj 


22 


3N- 


3x 


hjdXjd*, 


(3.14) 


Q 

This  results  in  a  six  by  six  element  stiffness  matrix,  S  ,  and 

e 

a  six  by  one  element  load  vector  F. 
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The  six  by  six  element  matrix  equations  of  equation  (3.14)  is 
wr 1 tten 


CSe]  Che  y  =  C F 6 > 


(3. 15) 


where 


1  = 


3  N  ,  3M.  9  N .  9  N  J 

11  «  .  22  „  .  \  a''  1  2’ 
9  K  i  9  *  i  9  x  2  9  x  2  j 


i,  j  =  1,2,3, .  .  .  6 


(3. 16) 


The  same  procedure  is  taken  for  all  elements.  Then  the  total 
stiffness  matrix  is  arrived  at  by  combining  the  element 
stiffness  matrices  using  structural  assembly  techniques 
(Zienkwicz,  1971) 


For  a  grid  of  elements,  a  system  of  banded  simultaneous 
equations  results: 


LSI  Chi  =  CF> 


(3. 17) 


where  CS1  is  the  global  matrix  of  cofficients  which 
incorporates  the  properties  of  the  materials  and  the  geometry  o+ 
the  elements;  Chi  is  the  vector  of  the  unknown  h’s  at  the 
nodes;  and  CFJ  is  the  load  vector.  After  the  introduction  of 
appropriate  boundary  conditions,  this  system  of  equations  can 


now  be  solved  by  computations  for  the  unknowns,  h,  yielding  the 
piezometric  head  distribution  in  the  solution  domain. 


-N  mV  *  V  .V  m  n  -  vW/JV  j»WV  *  .  *  - 


Free  Surface  Computation 


The  flow  domain  must  initially  be  established  assuming  an 
arbitrary  free  surface,  i.e.  the  user  initially  estimates  the 
free  surface  loacation.  All  boundary  conditions  are  solved  for 
the  preset  conditions  of  the  problem.  The  nodes  which  are  not 
fixed  are  allowed  to  move,  where  during  the  adjustment  the 
interior  nodes  are  also  moved  in  order  not  to  have  distortion 
of  the  elements.  Then  additional  iterations  are  made  until  all 
boundary  conditions  and  user  prescribed  error  correlations  have 
been  met  for  the  conf i gur at i on . 
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3.3  Variables  and  Dimensional  Analysis 

The  variables  investigated  in  this  study  are  shown  in  Figure 
3.1.  For  the  analysis  of  seepage  using  the  -finite  element 
program,  variables  were  changed  for  each  computer  run.  From 
each  run  data  was  generated  to  compare  the  changes  occur ing  for 
each  variable  change. 

Only  the  two-dimensional  flow  condition  was  considered  in  this 
study.  For  this  flow  condition,  the  building  drainage  system 
is  in  one  cr oss-sect i on  of  the  building  where  each  drain 
collects  the  seepage  and  collector  pipes  or  drains  take  the 
drainage  by  gravity  to  a  sump  or  some  exit  drain.  This  is 
depicted  by  Figure  3.4. 


I 


r.' 


The  total  discharge  of  the  system,  qT  ,  can  be  express  as  the 
following  function: 


—  ^(  Hr,  h  ,  dw,  S  p,  ® 


(3. 18) 


where 


q^  =  total  discharge  of  all  drains,  i.e.  the  sum  of 
discharges  through  each  seepage  face  for  the 
model,  ft  ^  /day/ft 

Hr  =  depth  of  saturated  so. 1  from  the  drain 
discharge  elevation  to  the  impervious 
boundary,  ft 

h  =  fixed  head  of  above  the  drains  ,  ft 
dw  =  drain  width,  ft 
Sp  =  drain  spacing,  ft 

K  =  hydraulic  conductivity  of  the  soil,  ft/day 
S  FS  =  difference  in  elevation  of  the  free  surface  to 
the  water  level  of  each  drain,  ft 
0  =  slope  of  the  impervious  boundary  with 

horizontal,  ft/ft 


Dimensional  analysis  of  equation  (3.18)  yields: 


H  d  S 


_  C  /  r  w  p 

-  o  ~~  1  r  >  ~  ^ 


h  h  h 


(3.19) 


Since  permeability  is  moved  to  the  left  side  of  equation  (3.19 
these  variable  is  made  proportional  to  the  discharge.  The 
elevation  of  the  free  surface,  Spg ,  is  an  output  variable,  so 


it  is  not  given  in  equation  <3.19.  Also  since  the  slope  of  the 
impervious  boundary,  0  ,  is  already  dimensionless  it  can  be 
given  as  a  single  term  in  FT/FT. 

The  effects  of  these  functional  variables  and  other  results  are 
given  in  the  form  of  dimensionless  plots  in  Chapter  4.  A 
discussion  of  the  results  is  given  in  that  chapter.  Since 
considerable  interest  is  in  quantities  of  discharge  and  the 
location  of  the  free  surface,  this  analysis  can  be  used  to  see 
trends  and  effects  that  each  variable  has  on  the  seepage 
domai n . 


3.4  Program  Verification 


To  verify  that  the  computer  program  used  in  this  study  is 
giving  valid  results,  two  other  methods  of  solutions  were  used 
for  a  selected  check  case.  Another  micro  computer  finite 
element  program,  SEEP,  and  an  appropriate  analytical  method 
based  on  Dupuit's  assumption,  was  used  for  the  simple 
conf iguration  of  a  drainage  blanket  with  a  single  pipe  type 
drain  (Duncan  and  Wong,  1  '85).  This  is  shown  in  Figure  3.5. 


Figure  3.5  Check  case  example  of  a  blanket  drain  with  single 
pipe  trench  drain. 


A  grid  was  generated  for  the  SEEP  program  to  match  the 
configuration  and  boundary  conditions  as  by  Figure  3.4.  The 
total  flow  for  this  check  case  is  calculated  by  the  SEEP 
program  as  based  on  the  Dupuit  assumption  of  a  fixed  domain  and 
fixed  amount  of  head.  Seep  also  tabulated  the  flow  quantities 
for  all  nodal  points  with  fixed  heads.  For  this  model,  the 


m 


mm 


■free  surface  converged  on  the  free  surface  location  with  an  error 


criteria  of  0.16  feet. 

The  location  of  the  free  surface  as  computed  by  the  Georgia 
Tech  model  converged  to  the  free  surface  having  an  error  factor 
of  only  0.001  feet,  a  much  more  accurate  and  sensitive 
convergence.  However,  the  two  models  compared  very  close  as 
far  as  the  location  of  the  free  surface.  Table  3.1  lists  the 
elevations  of  the  free  surface  nodes  for  the  two  models. 


TABLE  3. 1 


Elevations  o-f  the  Free  Sur-f  ace  Nodes 


SEEP: VP I 
Model 


Ga  Tech 
Model 


y,  coordinate 


y  coordinate 


0 

78.00 

78.00 

40 

72.55 

72.73 

70 

67.07 

67.57 

SO 

64.  27 

65.  17 

85 

62.03 

63.84 

87.5 

60.  39 

58.73 

90 

58 . 00 

58.  00 

92 

58 . 00 

58.  00 

95 

60.08 

58.32 

97.5 

61 . 16 

62.15 

1 00 

61.92 

62.52 

106 

63.  00 

63.28 

1  15 

63.  76 

63.93 

130 

64.  26 

64.42 

150 

64.42 

64.62 

drain  location 
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To  verify  that  -flaw  quantities  were  accurate  -for  the  Ga  Tech 
model,  results  were  also  compared  with  the  results  obtained 


from  a  Dupuit  assumption  solution,  developed  by  the  Federal 
Highway  Administration’s,  Highway  Subdrainage  Design  Manual , 
<1980).  Figure  3.6  was  used  to  calculate  the  quantity  of  flow 
to  the  drain.  The  following  calculations  based  on  Figure  3.6 
are  presented: 


b  1  L  .  91’ 

-  =  -  =  0.017  ,  -  =  -  =  1.57 

H  58’  H  58  ’ 

o  o 

K (H  -  H0)  0. 2(78’ -58’ ) 

-  = - =  2.17  (from  Figure  3.5) 

q  2  92 

0.2(20) 

therefore  q£  =  -  =  1-84  CFD 

2.  17 


q  ?  1 . 84 

K  0 .  2 


rp  '■n  <~\ 


(H  -  H  )2 

_ Q _ 

2 ( L  -b) 


(20) 


2(91-1) 


q  d 

—  =  9.22  +  2.22  =  11.44 

K 


therfore  qd=  0.2(11.44)  =  2.29  CFD 

The  results  of  this  computation  is  compared  with  the  results 
of  the  two  finite  elements  in  Table  3.2 


KV< 


Table  3.2  Discharge  to  the  drain  check  cases,  CFD 

Discharge,  CFD 

SEEP:  VF' I  Model  2.15 

Ga  Tech  Model  2.34 

FHA  Nomograph  2.29 

In  summary,  the  Ga  Tech  model  compares  very  well  with  other 
model  and  with  the  FHA  nomograph  -for  -flow  quantities.  The  Ga 
Tech  model  was  -found  to  more  error  criteria  sensitive.  The 
•free  surface  location  check  also  proved  to  compare  very  well 
with  the  SEEP  model. 
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CHAPTER  4 


RESULTS  AND  DISCUSSION 

4.1  Introduction 

This  chapter  discusses  the  analysis  and  results  of  the  two 
dimensional,  steady  state,  -free  surface  finite  element  study  of 
seepage  to  building  underdrain  systems.  Details  presented 
include  model  set-up  such  as  boudary  conditions,  grid 
development  and  basic  first  assumptions.  Output  variables 
considered  are  total  discharge  and  discharge  at  each  drain, 
plus  the  free  surface  elevation  between  the  drains. 

In  studying  sensitivity  of  the  model  for  the  input  variables, 
it  was  found  that  the  distance  from  the  Dirichlet  boundary 
face,  or  constant  head  boundary,  Li,  to  the  first  drain  was 
critical.  Li  was  found  to  be  critical  in  terms  of  output 
charcter i sti cs,  i.e.  discharge  and  location  of  the  free  surface 
between  drains.  Therfore,  two  sets  of  results  are  presented  for 
two  Li  distances  of  400  feet  and  80  feet. 

4.2  Model  Set-up 

Figure  4.1  shows  the  basic  grid  that  was  developed  for  this 
seepage  problem.  Using  150  elements  and  521  nodes,  the 
geometry  of  the  mesh  was  based  on  a  five  drain  system.  The 


sides  and  bottoms  of  the  drains  were  modeled  as  seepage  -faces. 


Using  the  two-dimensional  model,  the  -fi-fth  drain  is  considered 
to  be  at  the  building  centerline.  Assuming  symmetry,  only 
one-half  of  the  building  was  analysed.  The  left  most 
boundary,  as  shown  in  Figure  4.1,  is  a  Dirichlet  boundary 
condition  which  has  the  constant  head  required  to  allow  for  the 
seepage  to  enter  the  drains. 

Listed  below  are  some  basic  assumptions  for  the  study  made  for 
this  seepage  problem: 

<1)  Steady  state  flow  conditions 

(2)  Two-dimensional  analysis  where  the  right  most 
boundary  is  the  building  centerline  and  by  symmetry  the  flows 
are  only  one-half  of  the  total  for  the  total  width  of  the 

bui 1 di ng . 

(3)  Free  surface  nodes  move  to  the  computed  phreatic 
level,  i.e.  piezometric  head  equals  atmospheric  pressure  at  the 
upper  most  nodes 

(4)  Partial  saturated  flow  is  not  considered 

(5)  Evaporation  is  not  considered 

<6)  Retardation  is  not  considered 

(7)  Infiltration  is  not  considered 

(8)  The  soil  is  homogeous  being  either  isotropic  or 
an l sotrop l c . 
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Variables  Used 
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For  this  study,  key  variables  were  used  in  the  -finite  element 
model  to  look  at  their  effects  on  the  flow  of  seepage  for  the 
building  underdrain  system.  Two  major  design  parameters 
necessary  to  model  any  seepage  problem  is  the  total  discharge 
of  the  drain  system  and  the  location  of  the  free  surface.  For 
this  study,  the  following  output  variables  are  used: 


q^  =  discharge  or  seepage  quantity  to  each  drain 
(cubic  feet  per  day  per  linear  foot) 
q t  =  total  discharge  of  all  drains  (same  units) 

SpS  -  height  of  free  surface  above  drain  elevation. 


« 


Discharge  quantities  are  in  units  of  flow  quantity 
unit  length  of  the  total  flow  continum.  Therfore, 
total  system  should  be  multiplied  be  the  length  of 
structure,  i.e.  in  the  third  dimension. 


in  time  per 
flow  for  a 
the 


Input  variables  used  in  this  problem  to  study  their  effects  on 
seepage  are  listed  below: 


Kl 


ik  ■« 
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Li  =  distance  from  the  left  most  constant  head 
boundary  to  the  first  drain  (feet) 

K  =  hydraulic  conductivity  (feet /day) 

Kv  -  hydraulic  conductivity  in  the  vertical 
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direction  (feet/day) 


Fh  =  hydraulic  conductivity  in  the  horizontal 
direction  (-feet/day) 

h  —  hydraulic  head  of  flow  continum  above  the 
drain  elevation  (feet) 

Hr  =  depth  of  aquifer  between  the  elevation  of  the 
drain  and  the  impervious  boundary  (feet) 

Sp  =  spacing  between  drains  (feet) 

Dw  =  drain  width  (feet) 


8  =  slope  of  the  impervious  boundary  ('/.) 


For  better  reference  to  the  notation  of  the  variables  refer  to 
Figure  3.1  in  Chapter  3. 

For  this  finite  element  problem  the  units  of  feet  and  days  wer 
li  s  e  lI  t  a  g  i  v  e  r  esults  for  the  a  u.  t  p  u  t  v  ariable  s  i.  n  c  u  b  i  c  f  e  e  t  p  e 
day  and  feet  per  day  for  velocity. 

H  v  dra  u  lie:  c  o  n  duct  i  v  i  t  y  or  perm  e  a  b  i  1  i  t  y  v  a  1  u  e  s  u  s  e  d  f  o  r  t  h  i  s 
study  were  0.0002,  0.02,  0.2,  and  4  feet  per  day.  A  solid's 
hydraulic  conductivity  having  a  value  of  0. 0002  feet  per  day  i 
equivalent  to  a  fine  grained  soil  such  as  a  cl  ay .  Hydraul  i.  < . 
conduct i v i ties  are  tabulated  in  Table  4.1  to  show  cross 
reference  to  soil  types  and  permeabilities  in  units  a 
c e n  t  i.  m e t. e r  s  per  s e c o n d  . 


T  y  p  e  s 


Table  4.1  Hydraul ic  Conductivity  Values  For  Various  Soil 


Fine 

grained  soil 

0 . 0002 

f  eet /day 

7  m  10 

-7 

cm/ sec 

F i  ne 

grained  soil 

0. 002 

f  eet /day 

7  m  1 0 

^cm/  sec 

S  i  1 1 

0. 02 

feet /day 

7:;  10 

-7 

cm/ sec 

sand , 

silt 

0 . 2 

f  eet /day 

7x10 

-7 

cm/sec 

c  1  e  a  n 

sand 

4 

f  eet /day 

1 . 4x10 

-7 

cm/ sec 

The  other  input  variables  used  were  typically  realistic  values 
to  meet  typical  -field  conditions.  Listed  below  are  some  of 
this  values,  where  for  each  computor  run  specific  variables 
were  changed  and  the  others  held  constant.  The  intent  was  to 
see  what  effects  that  each  variable  change  made  on  the  model. 


Li  =  80’ 

,  400 

’  and  1000’ 

h  =10’ 

,  20  ’ 

and  40’ 

Hr  =  20’ 

,  58’ 

and  100’ 

Sp  =  15’ 

,  25’ 

and  35’ 

Dw  =  1 ’ , 

2’  and  4’ 

CD 

II 

rr  •/ 

U  /■  <i 

107.  and  207. 

Kv  =  .0002,  0.2,  4  and  15 


4.4  Model  Sensitivity 

The  model  was  developed  assuming  a  homogeneous  isotropic  media, 
except  tor  a  tew  cases  where  anisotropic  conditions  were  studied. 

The  nodes  ot  the  top  three  rows  ot  elements  were  allowed  to  move 
for  iteration  ot  the  tree  surface.  In  all  runs,  the  maximum 

I 

number  of  iterations  was  si;-:  and  on  an  average  three  to  four  were 
found  to  be  adequate.  The  model  was  made  to  iterate  to  an 
accuracy  of  0.001  feet  during  the  free  surface  moves. 

In  calibration  of  free  surface  models  careful  attention  must  be 
made  to  the  moving  nodes.  Problems  exist  when  the  elements  became 
less  than  zero  in  area  or  when  the  elements  actual  may  become  so 
distorted  that  they  may  become  inverted.  This  was  not  the  case 
for  the  mesh  used  in  this  study. 

In  studying  the  flow  quantities  and  velocities  for  the  basic 
geometries  and  variable  changes,  all  results  seemed  reasonable  and 
observed  trends  appeared  physically  possible.  However  large 
variation  of  output  was  noticed  by  the  change  in  distance  from  the 
constant  head  boundary  to  the  first  drain,  Li .  The  first  set  of 
results  is  presented  based  on  the  influence  length,  Li ,  being  10 
times  the  maximum  head  difference,  i.e.  10  x  40'  -  400’ .  Depth  of 

aqui fer  for  this  case  was  held  constant  at  100  feet.  Use  of  K' 
times  the  head  differential  is  based  on  Bear ’ s  (1979) 


recommendation. 


Figure  4.2  depicts  the  moved  -free  surface  corner  nodes  for  an 


£ 


$ 


I 


V\ 


« 

*> 


S 


s 

/J 


influence  length  of  400  feet.  For  comparison.  Figure  4.3  shows  the 
moved  free  surface  corner  nodes  for  Li  equal  to  00  feet  and  a 
thickness  of  aquifer.  Hr,  equal  to  58  feet. 
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are  present  ?d  to  show  the  dramatic  influence 


•  « 


S 


& 


I 


Figure  4.4  and  4.5 
that  the  distance  from  the  constant  head  boundary  from  the  first 
drain,  Li ,  has  on  total  discharge  and  the  maximum  free  surface 
elevation,  respectively.  For  further  reference,  the  total 
discharge,  q  ,  is  the  sum  of  the  drain  discharges  and  the  maximum 
free  surface  is  the  highest  elevation  of  water  above  the  drain 
elevation.  Lengths  between  300’  and  400’  appear  to  result  in  less 
change  in  flow  quantity  and  free  surface  changes.  This  equates  to 
the  distance  where  the  constant  head  boundary  away  from  the  drains 
begins  to  cause  the  least  effect  on  discharge  and  free  surface 
changes  is  appr ox i matel y  15  to  20  times  the  head  differential . 

This  perhaps  is  good  evidence  that  the  actual  field  conditons  far- 
each  particular  application  of  groundwater  movement  should 
carefully  be  looked  at  in  terms  of  this  distance. 
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4.5  Results 


This  section  will  present  the  computer  results  in  the  -form  of 
dimensionless  plots  and  other  nomographs  to  show  the  effects 
that  each  variable  change  has  on  the  seepage  problem.  Two 
sets  of  plots  are  given  where  the  first  set  is  for  the  constant 
head  boundary  distance,  Li,  equals  400  feet  and  the  depth  of 
rock,  Hr,  equals  100  feet.  The  second  set  of  plots  are 
essentially  the  same  as  the  first  except  Li  equals  80  feet  and 
the  Hr  equals  58  feet.  For  variable  notation  refer  to  Figure 
3.  1 . 


Drain  Discharge  (Li  =  400''  and  Hr  =  lOO*  ) 


Figure  4.6  shows  the  influence  of  drain  width,  Dw,  for  head 
d l f f er ent i al s,  h,  of  10,  20  and  40  feet  on  total  discharge  of 
the  drains.  This  shows  that  drain  width  is  not  particularly 
critical  for  low  head  differences.  It  is  reasonable  to  assume 
that  total  discharge  could  approach  a  constant  for  drain  widths 
greater  that  four  feet  and  head  less  than  20  feet.  Although 
total  discharge  increases  with  drain  width,  this  plot  does  not 
show  what  drain  width  is  the  most  effective  or  most  efficient. 
However,  the  influence  of  head  can  be  easily  observed 
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Figure  4.7  depicts  the  results  a-f  varying  anisotropic  soils 


with  varying  heads.  In  this  analysis,  vertical  hydraulic 
conductivity  varies  while  horizontal  conductivity  was  held 
constant.  Figure  4.7  shows  that  with  increase  in  vertical 
hydraulic  conductivity  discharge  increases.  This  appears  to 
satisfy  theory  ,  that  with  less  resistance  -flow  increases. 


Kv 

Kh 


Figure  4.7  Total  drain  discharge  vs.  anisotropy  with  varying 
head  differential. 

Li  =  400’,  Hr  =  100’,  Sp  =  15’,  Dw  =  2’ ,  Kh  =  0.2  ft/day 


Figures  4.B  through  4.12  are  graphs  which  show  the  discharge 
quantities  -for  each  at  the  tive  drains.  Figure  4.3  depicts  the 
effects  of  depth  of  impervious  layer,  Hr,  on  the  drain 
discharge  for  each  of  the  drains.  This  data  shows  that  with 
increased  depth  discharge  increases.  It  is  interesting  to  note 
that  practically  no  flow  is  recorded  for  drain  5  when  a  shallot- 
depth  of  impervious  layer  exists,  i.e.  Hr  =  20  feet. 
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Figure  4.3  Drain  discharge  vs.  thickness  of  aquifer. 
Li  =  400 ’ ,  h  =  20 ’ ,  Sp  =  25’,  Dw  =  2’,  K  =  0.2  ft /day 
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Figure  4.9  shows  the  effects  o-f  head  d  1  f  f  erent  i  al  on  the 
discharge  o-f  each  of  the  drains.  It  is  interesting  to  see  that 
though  the  first  drain  is  the  most  critical  drain,  i.  e.  a 
large  amount  of  the  total  flow  is  collected  in  this  drain,  the 
remaining  drains  collect  discharge  decreasing  linearly  from 
the  left.  This  is  possiblly  due  to  the  influence  of  the 
constant  head  boundary  distance,  Li,  being  400  feet  away.  This 
shows  that  for  the  drains,  flow  quantity  increases  linearly  with 
increase  in  head. 
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F i gure  4 . 9 
Li  =  400 ’ , 


Drain  discharge  vs. 
Hr  =  100’,  Sp  =  25', 
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Figure  4.10  depicts  the  results  o+  comparing  the  drain 
discharge  versus  the  spacing  between  drains.  It  is  o-f 
interest  to  see  that  drain  number  one  -for  spacing  equal  to  25 
-feet,  collects  more  seepage  than  the  other  discharges  -for  the 
other  spacings.  However,  at  drains  3  through  5  this  changes 
and  the  discharge  -for  15  Feet  spacing  is  greatest.  This  will 
be  analyzed  later  in  the  dimensional  analysis. 
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Figure  4.10  Drain  discharge  vs.  drain  spacing. 
Li  =  400’,  Hr  =  100’,  h  =  20%  Dw  =  2%  K  =  0% 
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Figure  4.11  shows  the  effects  of  anisotropy  of  the  soil 
matrix  on  drain  discharge.  It  is  of  interest  to  see  that 
the  discharge  for  Kv/Kh  equal  one  is  much  greater  than  the 
other  drain  discharges  for  Kv/Kh  less  than  one. 
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Figure  4.11  Drain  discharge  vs.  anisotropy. 

Li  =  400 7 ,  Hr  =  100',  Sp  =  15',  Dw  =  2,  h  =  20’ 
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Figure  4.12  depicts  data  generated  -for  the  condition  of  a 


sloping  impervious  layer.  The  slope  of  the  boundary  is  given 
as  a  percentage  of  fry-/ Ay.  The  depth  of  the  1  ayer  is  held 
constant  at  drain  one  and  the  1 ayer  slopes  clockwise  from  the 
positive  horizontal  axis. 

Though  actual  quantities  of  discharge  are  only  changing  in  very 
small  ammounts,  the  trend  is  seen  that  for  slopes  greater  than 
57.  the  discharge  increases  in  drains  3  through  5.  However  at 
some  slope  between  57.  and  107.  the  discharge  decreases.  It  is 
also  seen  that  the  discharge  of  drain  five  is  the  least  effected 
by  any  slope  change,  in  terms  of  discharge  quantity. 


1 


Free  Surface  Location  (Li  =  400 J  and  Hr  =  100J_) 

Figure  4-13  shows  the  influence  of  the  drain  width  on  the 

maximum  free  surface  height,  S  pc  -  In  all  cases  this  occured 

ramax 

between  drains  one  and  two.  This  plot  shows  how  the 
relationship  is  linear  and  that  the  height  of  the  free  surface 


is  not  greatly  affected  by  the  drain  width. 


In  Figures  4.16  through  4.20,  the  free  surface  corner  nodes  are 
plotted  for  the  different  cases  of  analysis.  To  see  the 


effects  of  the  free  surface  changes  versus  the  variable  changes 
only  nodes  between  the  drains  are  plotted.  In  the  sence  of 
design,  the  maximum  height  of  the  free  surface  would  be  of  the 
most  concern,  which  is  in  between  the  first  two  drains. 

However ,  it  is  of  interest  to  see  the  total  free  surf  ace  from  a 
standpoint  of  total  cross  section  of  the  building. 


In  Figure  4.16  the  effect  of  thickness  of  aquifer  on  location 
of  the  free  surface  is  shown.  It  is  of  interest  to  see,  as  in 
the  case  for  discharge  shown  by  Figure  4.8,  that  for  shallow  depths 
of  impervious  layer.  Hr  =  20’ ,  that  very  little  rise  in  the 
free  surface  occurs  in  between  drains  four  and  five. 


Figure  4.17  shows  the  effects  of  head  differential  on  the 
location  of  the  free  surface  nodes.  It  is  of  interest  to  see 
for  cases  of  head  equal  to  20  and  10  feet  the  free  surface 
nodes  are  approaching  gentle  mounds  of  water  in  between  drains. 
However  for  the  head  equal  to  40  feet  the  free  surface  peaks 
rather  abruptly  between  the  first  two  drains  but  resumes  the 
gentle  fall  in  free  surface  between  the  other  drains. 


Figure  4.18  shows  the  effects  of  drain  spacing  on  the  location 
of  the  free  surface  nodes.  This  shows  as  spacing  is  increased 
the  free  surface  increases. 


For  Figure  4.19,  the  e-f fects  o-f  the  sloping  aquifer  on  location 
of  the  free  surface  nodes  is  shown.  The  same  trend  is  seen 
here  as  in  discharge,  shown  by  Figure  4.12,  that  for  slopes 
greater  than  1 07.  the  free  surface  location  drops.  The  change 
in  free  surface  nodes  is  very  small  for  the  cases  of  slope 
equal  to  57.  and  107.. 


Figure  4.16  Free  Surface  corner  nodes  at  drains  with  effects 
impervious  boundary  layer. 

Li  =  400',  h  =  20',  Dw  =  2',  K  =  0.2  ft/day 


Surface  corner  nodes  with  drains  with  effects  of 
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Dimensionsal  Analysis  (Li  =  400“  and  Hr  =  100'') 

Figures  4.21  through  4.24  are  plots  of  the  data  generated  by  a 
dimensional  analysis  as  discussed  in  Chapter  3.  The  discharge 
values  given  are  unit  discharge  quantities  and  each  case  is  -for 
the  particular  variables  noted  in  each  -figure.  Unit  discharge 
is  defined  as  the  sum  of  the  total  drain  discharges  for  each 
of  the  drains  per  unit  width  in  the  z  axis  or  the  third 
di  men si  on . 


In  Figure  4.21  the  effect  of  drain  width  is  analyzed.  This 
plot  shows  that  with  increased  drain  width  discharge  increases, 
but  with  increased  head  the  discharge  may  increase  at  a  much 
faster  rate.  From  this  analysis  it  would  appear  that  the 
greatest  increases  of  flow  rate  would  be  for  the  drain  widths 
between  one  and  two  feet  wide.  The  following  analysis  is 
perf armed: 


l ncrease 


Drain  Widths 

head  123  4 

10  0.034  0.053  0.065  0.073 

—  0.019  0 . 0 12  0 .  009 

20  0.035  0.055  0.070  0.086 

increase  - 0.02  0.015  0.016 

4 0  0 . 0 3 6  0 .057  0.075  0 . 0 9 

increase  - *-  0.021  0.018  0.015 
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This  analysis  shows  the  increases  of  the  -flow  -function  -for  the 
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di-f-ferent  head  increases  as 
seen  -for  the  values  between 
Therfore  the  most  effective 
widths  between  1  and  2  feet, 
stated,  that  for  Dw/h  ratios 
increases  at  a  rapid  rate. 


shown.  The  maximum  increase  is 
drain  widths  of  one  and  two. 
increase  in  flow  occurs  for  drain 
Another  relationship  could  be 
less  than  0.2  the  flow  quantity 
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In  Figure  4.22  the  effects  of  thickness  of  aquifer  is  shown  for 
the  dimensional  analysis.  This  plot  concludes  that  for  Hr /h 
ratios  greater  than  2  or  3  the  net  increase  in  flow  is 
approaching  a  linear  rel ati onshi p.  This  states  that  with 
increasing  thickness  of  the  aquifer  the  dicharge  is  expected  to 
increase  very  near  a  linear  increase  and  for  very  shallow 
thicknesses  the  flow  becomes  very  small. 


Figure  4.22  Dimensional  Analysis,  Depths  of  impervious  layer. 
Li  =  400 ’ ,  Sp  =  25’,  Dw  =  2’  ,  K  =  0.2  ft/day 
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The  effects  of  drain  spacing  is  presented  in  Figure  4.23.  The 
•flow  function  for  this  analysis  is  treated  differently  for  the 


spacing  changes.  As  the  spacing  was  increased  for  each  of  the 
computer  runs  the  flow  domain  increased  accordingly  to 
accommodate  the  five  drains  at  the  appropriate  spacing. 
Therfore  to  gain  understanding  of  the  effects  of  flow 
relationships  with  respect  to  drain  spacing,  each  term  was 
normalised  by  the  cor respond i ng  spacing  between  the  drains  for 
each  head  increase.  This  analysis  shows  for  decreased  spacing 


discharge  increases  and  as  head  increases  so  the  discharge. 


Figure  4.23  Dimensional  Analysis,  Drain  spacing. 

Li  =  400,  Hr  =  100’  ,  Dw  =  2’ ,  K  =  0.2  ft/day 

Figure  4.24  shows  the  dimensional  analysis  for  effects  of  the 
sloping  aquifer.  For  this  particular  flow  matrix  as  modeled, 
a  definate  break  in  flow  increase  is  noticed.  The  peak- 
discharge  function  i s  at  slope,  0  equal  0.05. 


Results  For  Li 


80'’  and  Hr 


58’ 


The  following  plots  and  nomographs  are  the  results  for  the 
computor  runs  where  distance  from  the  constant  head  boundary  to 
the  first  drain,  Li,  equals  80  feet  and  the  thickness  of  the 
aquifer  equals  58  feet.  The  plots  are  essentially  the  same  as 
for  the  previous  case  except  these  variables  are  held 
constant . 


Drain  Discharge  <  Li  =  80’  and  Hr  =  58’  ) 


Figure  4.25  shows  the  effects  of  drain  width  on  total  drain 
discharge.  The  results  are  generally  the  same  as  shown  in 
Figure  4.6  except  discharge  quantities  range  from  3  to  4  times 
greater . 
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The  results  of  anisotropy  of  the  soil  continum  is  shown  by 


Figure  4.26  for  the  short  Li  distance  of  80  feet.  The  trend  is 
slightly  different  than  for  the  case  for  Li  equal  400  feet  as 
shown  by  Figure  4.7.  For  the  largest  head  value  of  40  feet  the 
discharge  decreases  just  slightly  with  decrease  in  vertical 
hydraulic  conducti vi ty,  but  for  a  Kv/Kh  ratio  less  than  0.1  the 
discharge  decrease  is  very  dramatic.  However  for  the  low  head 
value  of  10  feet  the  general  decreasing  trend  is  prevelant  as 
the  case  for  Figure  4.7. 

For  this  case  the  flow  quantities  have  increased  approximately 
4  to  5  times.  This  shows  the  large  effects  that  the  distance 
of  the  constant  head  boundary  has  on  the  seepage  quantities. 
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Figures  4.27  and  4.28  depict  the  same  general  trend  o-f  maximum 

| 

•flow  at  drain  one  and  a  general  decrease  in  -flow  -for  each 

*  consectutive  drain  is  noticed  -for  effects  of  thickness  of 

* 

* 

aquifier  and  head  respectively.  Flow  quantities  are  again  3  to 
!  4  times  greater  for  this  case  than  for  the  earlier  case  as 

shown  by  Figures  4.8  and  4.9.  For  this  case  of  shorter 
distance  to  the  constant  head  boundary  equal  80  feet,  the 
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Figure  4.27  Drain  discharge  vs.  depth  of  aquifer 
Li  =  80’  ,  h  =  20’ ,  Sp  =  25’,  K  =  0.2  ft/day,  Dw  = 


The  effects  of  drain  spacing  for  distance  of  the  constant  head 


m 


■ 

SPACING 

15' 

1 

SPACING 

25' 

a 

SPACING 

35' 

In  Figure  4.30  the  effects  of  slope  o-f  the  impervious  boundary. 


y  ,  is  shown  tor  drain  discharge.  A  significant  change  is 
noticed  for  this  plot  when  compared  to  the  earlier  plot  in 
Figure  4.12  for  sloping  aquifer.  For  Li  equal  to  400  feet,  the 
peak  discharge  occured  for  a  slope  between  5  and  10  percent. 

In  this  case,  where  Li  equals  80  feet,  the  peak  flow  is  for  the 
case  of  the  steepest  slope  equal  to  20  percent.  For  drain  one, 
the  flow  is  8  times  greater  than  the  flow  for  drain  one  as 
shown  by  Figure  4.12,  Li  =  400  feet.  The  flow  is  5  times 
greater  in  drain  five  in  this  case  than  the  flow  in  drain  five 
for  the  case  of  Li  =  400  feet. 
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Free  Surface  Location  (  Li  =  SO’  and  Hr  =  58 '  ) 

Figures  4.31  through  4.36  show  the  results  o-f  the  -free  surface 
location  plots  -for  the  case  o-f  Li  equal  30  -feet  and  Hr  equal  58 
feet.  The  same  general  trends  are  shown  for  this  plots  as 
shown  by  the  earlier  case  of  Li  equal  400  feet  and  Hr  equal  100 
feet.  In  general ,  the  free  surface  has  moved  up  approximately 
three  times  in  height.  An  exception  to  the  earlier  case  of  is 
shown  in  Figure  4.36.  This  figure  shows  the  location  of  the 
free  surface  nodes  with  the  effects  of  slope  of  the  impervious 
boundary.  It  is  somewhat  different  in  this  case  than  in  the 
earlier  case  as  shown  in  Figure  4.19  for  Li  equal  400  feet.  In 
this  case  the  free  surface  continues  to  increase  in  height  as 
the  slope  of  the  impervious  boundary  increases. 
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Figure  4.32  Maximum  free  surface  height  v^.  -train  spacing 
ith  head  varying. 

i  =  80’ ,  Hr  =  58’,  Dw  =  2’ ,  K  =  0.2  ft/day 
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Free  surface  corner  nodes  with  effects  of 
Depth  of  impervious  boundary. 

Li  =  80' ,  h  =  20' ,  Dw  =  2' ,  K  =  0.2  ft/day 
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Figure  4.36  Free  surface  corner  nodes  with  effects  of 
slope  of  impervious  boundary. 

Li  =  80',  Hr  =58' (first  drain),  Dw  =  2',  h  =  20',  K  =  0.2  ft/day 


Dimensional  Analysis  (  Li  =  30’  and  Hr  =  58'  ) 

The  e-f  +  ects  oT  drain  width  as  shown  by  the  dimensional  analysis 
is  shown  by  Figure  4.37.  The  general  trend  o-f  the  discharge 
increasing  with  drain  width  is  seen  in  the  data.  When  compared 
to  the  earlier  case  -for  Li  equal  to  400  -feet,  it  is  seen  that 
total  discharge  is  greatly  e-f-fected,  by  a  -factor  o-f  three 


_Si 

Kh 


Dw/h 

Figure  4.37  Dimensional  Analysis,  drain  width. 
Li  =  80',  Hr  =  58',  Sp  =  25',  K  =  0.2  ft/day 


Figure  4.38  shows  the  dimensional  analysis  results  -for 
thichness  o-f  the  aquifer.  In  this  case,  Li  equal  60  feet,  the 
relationship  is  more  linear  than  in  the  case  of  Li  equal  400 
feet  as  shown  in  Figure  4.22.  This  results  in  the  conclusion 
that  the  effects  of  Hr  has  a  pronounce  influence  on  the  flow 
quantities  when  coupled  with  the  distance  Li.  For  a  shorter 
distance  of  Li,  the  effects  of  depths  of  Hr  is  more  crictical 
and  more  of  a  linear  relationship.  However  when  the  distance 
from  the  constant  head  boundary  is  larger,  as  in  the  earlier- 
case,  the  effects  of  Hr  are  not  as  effective  when  Hr  is  very 
smal 1  . 
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Dimensional  Analysis,  depth  o-f  aqui-fer 
=  25%  Dw  =  2%  K  =  0.2  ft /day 


Figure  4.39  shows  the  results  of  the  dimensional  analysis  for 
drain  spacing.  Again  the  ordinate  is  normalised  with  respect 
to  the  area  subtended  by  each  flow  continum  for  each  case  of 
spacing.  Flow  quantity  would  equate  to  an  increase  of 
approximately  5  to  6  times  from  this  case  of  Li  equal  80  feet 
versus  Li  equal  400  as  shown  in  Figure  4.23. 
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Figure  4.39  Dimensional  Analysis,  drain  spacing, 
i  i  .  Hr  =  50'.  Dw  =  2’.  K  =  0.2  ft/day 


4.<b  Summary 


This  chapter  has  included  the  results  o-f  many  computer  runs  -for 
the  analysis  of  a  building  underdrain  system.  The  importance 
o-f  model  sensitivity  must  be  stressed,  in  that  key  variables  may 
effect  the  results  very  significantly.  In  this  model,  it  was 
found  that  the  most  crictical  variables  affecting  the  total 
discharges  of  the  drains  and  location  of  the  free  surface  was 
the  distance  from  the  constant  head  boundary  to  the  drain 
system,  Li,  and  thickness  of  the  aquifer,  Hr. 

The  results  obtained  appear  to  be  reasonable  and  within  limits 
of  possible  field  occurances.  However,  using  these  results  for 
quantifying  flows  and  the  free  surface  locations  requires  a 
very  careful  study  and  background  investigation  of  all 
geometric  conf i gur at i ons  and  soil  conditions.  It  is  intended 
that  this  study  has  presented  the  results  in  a  manner  that  one 
can  observe  the  effects  of  the  variables  studied  on  such  a 
underdrain  system. 

In  summary,  the  following  observations  are  discussed  for  the 
variables  used  in  this  analysis. 

Constant  Head  Boundary  Distance,  Li 

As  distance  from  the  constant  head  boundary  to  the  system 


increases  the  flow  quantities  increase  -for  the  drain  system. 

In  terms  of  a  correlation  with  hydraulic  head  differential 
above  the  drain,  as  used  by  other  authors,  this  model  study 
shows  that  for  a  distance  of  approximately  15  times  the  head, 
discharge  appears  to  decrease.  In  terms  of  gradient,  this 
equates  to  a  1  over  15  ratio  or  0.067  gradient.  A  slightly 
higher  distance  was  observed  in  the  analysis  of  this  influence 
to  location  of  the  free  surface. 

Depth  of  Aquifer,  Hr 

The  results  show  that  with  increased  aquifer  thickness  the 
drain  discharge  increases,  increasing  very  near  to  a  linear 
rel ati onshi p.  It  was  also  observed  that  for  the  very  shallow 
thicknesses,  the  last  drains  collected  very  little  flow.  The 
height  of  the  free  surface  is  affected  in  the  same  fashion. 

With  increased  thickness  of  aquifer  the  height  of  the  free 
surface  increases  between  drains. 

Drain  Spacing,  Sp 

Drain  spacing  effects  on  discharge  was  greatly  affected  by  the 
distance  of  the  system  from  the  constant  head  boundary.  A 
substantial  increase  of  5  to  6  times  the  discharge  was  observed 
for  the  two  distances  used  in  the  results  presented.  It  was 
observed  that  as  head  increases,  discharge  increases  for  all 
drain  spacings;  and  as  drain  spacing  increases,  discharge 
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Drain  Width,  Dw 

With  increased  drain  width  the  discharge  increases  -for  all 
other  variables  constant.  However,  there  must  be  a  trade  o-f-f 
-for  e-f-ficiency  o-f  the  complete  system.  It  appears  that  a  drain 
width  between  1  to  2  feet  is  the  most  effective  width.  This  is 
also  a  plus  from  a  standpoint  that  lateral  seepage  can  be 
easily  drained  into  thinner  trenches  or  deeper  constructed 
drains. 

Slope  of  Impervious  Boundary,  0 


For  the  distance  of  the  constant  head  boundary  to  the  drain 
system  equal  to  80  feet,  the  drain  discharge  was  observed  to 
increase  with  increased  slope.  However,  for  the  distance,  Li 
equal  to  400  feet,  the  maximum  discharge  occurs  at  a  slope 
between  5  to  10  percent.  This  occurance  is  not  fully 
understood  and  perhaps  should  be  further  studied  to  develope 
some  reasoning.  This  requires  that  careful  consideration  of 
specific  boundary  geometries  be  used  in  model  studies  for  seepage 
anal ysi s. 
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CHAPTER  5 


CONCLUSS I DN 

In  this  day  and  age,  more  and  more  sites  are  chosen  tor 
buildings  and  engineered  structures  that  are  less  desirable  to 
work  with  from  a  standpoint  of  site  conditions,  poor  soils  and 
water  problems.  A  less  desirable  site  may  easily  be  where  the 
groundwater  table  is  very  near  to  the  surface.  Thertore  the 
building’s  design  must  incorporate  some  thought  into  how  the 
groundwater  will  be  handled  and  what  effects  it  may  have  on  the 
str ucture. 

Engineers  faced  with  the  analysis  of  seepage  for  buildings  to 
be  built  below  the  groundwater  table  must  consider  such 
conditions  of  uplift  pressures,  critical  velocities,  seepage 
control  of  the  water  itself  or  particle  migration.  Hopefully 
the  design  will  enhance  stability  of  the  structure  and  prolong 
the  usefulness  and  life  of  the  structure. 

In  performing  the  literature  review  on  this  topic,  it  was  found 

that  many  authors  and  researchers  are  more  and  more  relying  on 

and  developing  computor  models  to  perform  the  analysis  of 

cample:;  seepage  problems.  With  the  use  of  general  seepage 

theory  and  mathematical  representation  of  the  seepage  theories, 

119 


i 


f 

v> 

.V 

,v 

■A. 

A. 

y, 

i 


4%‘ 

.V 


! 


,  J 
9J 


I 


tarr  <V  *,  '  -  « 


the  computer  programs  can  quickly  and  efficiently  aid  the 
endgineer  in  studying  the  flow  conditions.  As  in  this  study, 
the  analysis  of  one  particular  seepage  problem  using  a  finite 
element  program  made  it  possible  to  study  many  key  variables 
and  develope  many  re  1  at i onsh 1 ps . 

For  this  study,  underdrain  systems  for  buildings  built  below 
the  groundwater  table  were  researched.  By  using  a  finite 
element  program,  r-nd  by  developing  a  mesh  for  this  particular 
flow  matrix,  many  variables  of  the  analysis  were  studied.  If 
was  found  that  the  model  could  easily  be  used  to  etn  ienM, 
and  accurately  study  these  variables.  ♦  ma  u  e-  .  *  •  e 

dicharge  of  the  svstem  and  1  oc  a  t  1  •  <  ♦  *  -  e  *  •  ee  •  •  <  .  ~  • 

studi  ed  in  detail  versus  f  f-e  .  <*«•  . 

This  r  esear  c  h  h  a  s  a  *  *  e  r.  [  *  :  • 

and  usef  u  1  cr  hem  ^ 
seepage  problem.  •  ■. 
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APPENDIX  A- 1 


USER’S  MANUAL  FOR 


ARAL  SEEPAGE  PROGRAM 


I.  INTRODUCTION: 


This  seepage  program  is  a  -finite  element  model  for  solving 
unconfined  or  confined  two  dimensional  or  axisymetric  seepage 
problems.  For  unconfined  flow  analysis  the  model  iterates  by 
moving  the  free  surface  nodes  to  the  phreatic  level. 


In  general,  the  program  allows  for  the  porous  soil  mass  to  be 
divided  into  discrete  elements  and  nodes  using  either  a  vertical 
plane  model  (two-dimensional  ,>:-y)  or  axisymmetric  grid  structure. 
Then  applying  known  boundary  conditions,  the  program 
simultaneously  solves  for  each  element,  the  continuity  equation 
which  satisfies  Darcy’s  Law. 


Figure  1  shows  a  very  simple  case  of  a  two  dimensional  unconfined 
steady  state  condition  which  can  serve  as  an  example.  In  this 
example  the  free  surface  location  is  unknown  and  the  user  must 
begin  the  process  by  i ni t i al 1 i z i ng  the  input  by  approx i mat l ng  the 
location.  Then  the  program  will  iterate  to  the  free  surface  level 
until  the  free  surface  converges  to  satisfy  the  equations  of 
f  1  ow . 


II.  PROGRAM  EXECUTION  AND  ANALYSIS: 


The  user  has  the  option  of  performing  a  two  dimensional  problem  or 
an  axisymmetric  problem.  The  two  dimensional  case  simply  show  b 
figure  1  will  serve  as  an  example  to  explain  the  parameters.  An 
axisymetric  application  might  be  a  well  source  problem  looking  at 
stratified  soils  having  varying  permeabilities  at  different 
elevations  through  the  porous  media. 

III.  BOUNDARY  CONDITIONS  OF  FLOW: 

The  potential,  0,  of  flow  must  be  defined  as  some  directional 
f unct i on 


0  =  ^  ( >: ,  y ,  z  ,  t ) 

A  Dirichlet  boundary  condition  is  where  the  piezometric  head,  h  is 
constant,  or 


h 


.y»z 


» 


t) 


is  solved  for  at  a  particular  point.  This  boundary  condition 
occurs  whenever  the  flow  domain  is  adjacent  to  an  open  body  of 
fluid.  Since  the  piezometric  head  is  constant  at  the  open  body  of 
water,  the  piezometric  head  of  the  porous  media  at  the  interface 


will  be  constant.  This  boundary  condition  denotes  a  non-  zero 
■flux  or  -flow  is  to  be  considered  possible. 

Another  boundary  condition  -for  -flow  potential  is  the  Neuman 
boundaries.  This  is  the  case  where  -flow  cannot  move  through  the 
porous  media  at  a  specific  point  defined,  i.e.  an  impervious  layer 
of  stratum. 

IV.  FINITE  ELEMENT  GRID  NOTATION: 

The  porous  soil  mass  to  be  analysed  must  be  divided  into  a 
grid  of  elements  and  nodes.  The  program  uses  any  shape  second-order 
i soparametr i c  element,  shown  by  figure  2  as  a  quadral ateral 
element  containing  four  corner  nodes  and  four  midpoint  nodes. 


Figure  2  Element  and  Node  Notation. 


The  nodes  must  be  numbered  -from  some  beginning  corner,  preferably 


the  number  one  node  should  be  at  the  bottom  left  corner  of  the 
mesh.  Then  node  numbering  increases  to  the  end  of  each  column, 
beginning  at  the  bottom  of  each  column  throughout  the  mesh.  One 
can  see  the  numbering  scheme  by  the  examples  presented.  Elements 
are  numbered  the  same  as  the  nodes. 


The  program  automat i cal  1 y  generates  nodes  at  equal  spacing  between 
corner  nodes  which  are  defined  by  the  user.  This  feature 
requires  that  careful  attention  be  made  to  the  node  numbering  scheme 
once  corner  nodes  have  been  established.  This  feature  is  good  for 
creating  finer  mesh  grids  in  particular  areas  of  interest. 

Base  elements  are  used  to  create  stationary  reference  elements 
typically  at  the  bottom  of  the  mesh  or  away  from  the  free  surface 
elements.  These  elements  are  also  used  to  define  changes  in 
hydraulic  conductivity  in  the  x  and  y  di recti ons (Kx  and  Ky). 


The  free  surface  or  phreatic  surface  of  a  media  is  defined  by  the 
user  as  a  first  approximation  of  it’s  actual  location.  In  the  case 
of  a  cutoff  wall,  shown  in  Figure  3,  there  is  no  free  surface  and 
the  seepage  model  solves  for  the  flow  character i st i cs  of  each 
element  and  node.  However,  in  the  case  of  a  dam  as  shown  in 
Figure  3,  the  free  surface  is  unknown  and  the  program  will  solve 

for  the  phreatic  surface  where  the  fluid  pressure  i s  at 

A- 6 


atmospheric  pressure. 


In  confined  aquifers  the  water  pressures 
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would  possible  be  higher  than  atmospheric  pressure. 
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Seepage  faces  are  defined  by  the  user  to  denote  a  boundary  or 
where  seepage  is. to  be  allowed.  Here  flow  is  qunatified  and  the 
model  solves  for  the  velocity  and  discharge  at  this  face. 


Figure  3  Example  of  confined  flow  and  unconfined  flow. 
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Permeability  o-f  the  porous  media  often  depends  on  the  direction 


of  flow.  For  the  condition  of  soil  matrix  being  anisotropic  the 
velocities  and  gradients  are  generally  not  parallel.  For  this 
condition  the  velocities  can  be  written  in  the  Darcy  Law  form: 


=  — K>:  >: 


c)h 

-Kyy  ___ 


=  -Kzz 


Normally  the  hydraulic  conductivity  of  teh  media  is  similar  in  the 
y  direction.  Therfore  treatment  of  this  normal  condition  requires 
use  of  only  two  direcetions.  For  this  simplification,  the  program 
uses  the  notation  of  the  two  dimensional  case.  Figure  4  shows  the 
notation  that  should  be  used  for  describing  the  hydraulic 
conductivity  conditions. 


Figure  4  Hydraulic  Conductivity  Notation. 
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INPUT  VARIABLES  AND  DATA  INPUT 
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Input  can  be  given  in  unformated  notation  where  each  data  input  is 
separated  by  a  comma  or  space.  Careful  attention,  however,  must 
be  made  to  use  the  apppropriate  number  input  data  per  line  as 
described  by  this  guide. 


Title  and  Type  Problem  Card. 

TITLE  The  title  can  be  up  to  72  characters. 

DIMEN  Flow  problem  e  problem  is  a  two  dimensional  analysis 
input  2  and  for  an  axisymetric  problem  input  1. 

Plane  Flow  Problem  2 

Axisymetric  Problem  1 

FERR  Having  performed  a  sensitivity  of  the  model  or  decided 
on  an  error  critiria,  input  the  error  required. 


Control  cards 


NNF'C  Total  number  of  corner  nodal  points. 

These  nodes  are  used  to  define  corners  of  elements 
which  are  defined  by  input  coordinates.  Nodes  that 
are  not  listed  will  be  generated  between  two  corner 
nodes  for  each  vertical  column  or  row  deliniation. 


NELEMC  Total  number  of  base  elements. 

Base  elements  are  typically  at  the  bottom  of  a  series 
of  elements  that  change  the  hydraulic  conductivity  of 
the  elements  above  these  base  elements.  These  elements 
are  established  as  reference  elements  for  the  program 
to  describe  permeability  of  each  element. 

NF'BQC  Total  number  of  Neuman  Boundary  Faces. 

A  Neuman  face  is  defined  as  a  line  of  nodes  that  may 
create  a  barrior  or  impervious  wall  for  the  flow.  If 
two  corner  nodes  define  a  line  of  several  nodes  that 
makes  an  impervious  barrier,  this  is  one  Neuman 
Boundary  face. 

NPFS  Total  number  of  corner  nodes  on  the  free  surface. 

The  user  must  specify  particular  corner  nodes  which 
reflect  the  first  approx i mat i on  of  the  unknown  free 
surface.  The  program  uses  these  nodes  and  the  midpoint 
nodes  to  converge  to  the  free  surface. 
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ITGIV 


C-f  r 


Number  of  iterations  desired. 

This  is  a  safety  feature  tor  possible  lupi.t 
computor  time  limitation. 
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NPBDC  Total  number  ot  Dirichlet  Boundary  tac.es. 

Flow  is  allowed  at  these  faces.  Do  not  count  the  m d 
point  nodes  as  part  of  the  face. 


NSPFAC  Total  number  of  seepage  faces  with  Dirichlet 

boundaries.  For  each  face  defined  by  two  corner  ru  des 
for  the  NPBDC  input  variable  a  NSPFAC  value  should  be 
given.  As  in 


NPBDC  =  2 
NPBDC  *  3 


NSPFAC  =  1 
NSPFAC  =  2 


ITIME 


Time  variable.  Input  0  =  time  independent. 

This  is  used  for  transient  computation  and  is  not 
used  in  steady  state  computations. 


Corner  Node  Location  Cards 

Corner  node  number  <n) 

(n>  >i  coordinate  of  node(n) 

y(n)  y  coordinate  of  node(n) 

Note:  There  should  be  NNPC  number  of  these  cards.  Use 

appropriate  dimensions,  i.e.  feet,  meters  which  will  be  carried 
throughout  the  remaining  input  variables. 

NPMIS  Enter  a  "1"  or  "0"  where 

1  -  additional  nodes  are  to  be  generated  after  this 
node. 

0  -  no  nodes  will  be  generated  after  this  node. 
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Base  Reference  Elements  Identification  Cards 
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NOD 


For  each  base  element (N),  give  a  counter-clockwise 
numbering  for  each  node  in  the  element. 
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NMIS  Total  number  of  elements  above  the  base  element  before 
the  next  base  element. 

XK(N)  Hydraulic  conductivity  of  the  elements  defined  by  the 
base  element  in  the  x  direction. 

YK(N)  Hydraulic  conductivity  in  the  y  direction. 

ALP (N)  Angle  of  possible  strati f i cati on  or  plane  of  soil 
matrix  with  the  horizontal. 


Note:  There  should  be  NELEMC  number  of  these  cards  provided. 

Neuman  Boundary  Condition  Cards 

NSTART  The  beginning  or  lowest  corner  node  number  of  each 
Neuman  Boundary  defined. 

BVAL  The  Neuman  Boundary  value  that  defines  the  known 
potential  height  or  water  level  measured  from  the 
datum  reference  line  of  the  coordinate  system. 

NBSAME  Total  number  of  boundary  faces  having  same  Neuman 

boundary  condition.  Each  face  is  a  line  made  by  two 
corner  nodes  of  the  same  Neuman  boundary. 
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NINC  The  number  of  nodal  number  increases  or  numerical 

difference  to  the  next  corner  node  from  node  NSTART. 
Do  not  count  mid-nodes  only  corner  nodes. 

Note:  There  should  be  NF'BOC  cards  provided  for  this  set. 


Dirichlet  Boundary  Condition  Cards 

NSTAR  The  beginning  or  lowest  corner  node  number  which 
defines  each  Dirichlet  boundary. 

DBVAL  The  Dirichlet  boundary  value  that  defines  the  potential 
or  water  level  height  measured  from  the  data  reference 
elevation  of  the  coordinate  system. 

NDSAME  Total  number  of  boundary  faces  of  this  Dirichlet 
boundary  condition.  Each  face  is  counted  between 
corner  nodes  defining  the  boundary  or  number  of 
elements. 

NDINC  The  number  of  nodal  number  increases  or  numerical 

difference  to  the  next  corner  node  from  node  NSTAR. 

NMID  The  number  of  nodal  number  increases  from  the  NSTAR  to 
the  next  mid-node. 

Note:  There  should  be  NPBDC  cards  provided  for  this  set. 


Free  Surface  Nodal  Cards  (Maximum  16  pairs  per  card) 

NF'FSA  Free  surface  corner  node  number. 

NF'BOT  Node  number  of  the  bottom  of  the  movable  section  of 
the  grid  that  defines  the  limits  of  moving  nodes  for 
NF'FSA. 

Note:  There  should  be  NF'FS  pairs  of  node  numbers  given  for  a 

maximum  of  16  pairs  per  card.  Use  as  many  cards  as  required.  The 
bottom  node  (NF'BOT)  defines  the  bottom  of  a  seepage  face 
correspondi ng  to  the  grid  geometry  and  corner  node  defining  the 
free  surface.  The  program  uses  this  band  of  nodes  and  elements  to 
converge  on  the  free  surface. 


Mid-point  Free  Surface  Nodal  Cards. 
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ft 

NSINC 

NF’FSM  The  mid-point  nodal  numbers  between  corner  nodes 
defining  the  free  surface  (NF'FS  -  1). 


Seepage  Face  Cards 

NSF'C  Beginning  node  number  of  the  defined  seepage  face. 


seepage  face  NSF'C. 

NSINC  Node  numbering  increase  between  nodes  defined  by 
NSSAME. 

Note:  There  should  be  NSF'FAC  cards  provided  for  this  set. 

Statement  and  End  Card 

End  of  Problem  Statement  should  be  less  than  72  characters  and 
end  with  a  period. 

ND I MEN  Input  a  "3"  to  end  the  problem. 

ERR  Specify  any  error  as  for  FERR. 
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UNCONFINED  FLOW  EXAMPLE  PROBLEM 
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Title  and  Type  Problem  Card 

TITLE  -  SILTY  SAND  DAM 

DIMEN  -  2  (two  dimensional) 

FERR  -  O.Ol  (error  criterion  is  0.01) 

Control  Cards 

NMPC  -  12  (corner  nodes  are  18,  52  35,  45,  52,  62,  69,  79,  86  and 
96) 

NELEMC-  5  (  five  base  elements  -  shaded) 

NPBOC  -  0  (  no  Neuman  Boundaries) 

NPFS  -6(6  nodal  points  on  the  free  surface  -  11,  28,  45,  62, 

79  and  96) 

ITGIV  -  50  (  maximum  number  of  iterations) 

NPBDC  -  2  (There  are  two  faces  with  Dirichlet  boundary  conditions, 
the  left  side  and  the  right  side  of  the  mesh  which 
forms  the  dam) 

NSF'FAC-  1  (There  is  one  seepage  face  -  on  the  right  side  where 
flow  is  to  be  considered.) 

ITIME  -  0  (Time  independent  for  steady  state  conditions) 
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Corner  Node  Location  Cards 

1,0, 0 ,  1 

Where,  the  -first  node  is  number  1  having  an  ::  coordinate  of  0  and 
a  y  coordinate  o-f  0.  Generation  of  additional  nodes  will  follow. 

1 1 , 0 ,  110,0 

Where,  node  11  is  at  the  end  of  the  first  column  at 
coordinates  0  and  110  and  no  generation  is  required  for  other 
nodes  in  this  column. 


18, 

200, 

0,  1 

28, 

200, 

99, 

o 

35, 

400, 

0,  1 

45, 

400, 

88, 

0 

52, 

600, 

0,  1 

62, 

600, 

97, 

o 

69, 

800, 

0,  1 

29, 

800, 

66, 

0 

86, 

1 000 

,0, 

1 

96, 

1 000 

cc 
^  JJ 

Base  Reference  Elements  Identification  Cards 
1,12, 18, 19, 20, 13, 3,2,4,0.002,0. 002, 0 

Counter  clockwise  numbering  of  nodes  for  each  element  and  the 
permeabi 1 i ty ,  0.002  ft/d ay  for  this  base  element  and  4  elements 
above  the  base  element. 

1 8 , 29 , 35 , 36 , 37 , 30 , 20 , 1 9 , 4 , 0 . 002 , 0 . 002 , 0 
35, 46,52,53,54, 47,37,36, 4,0.002,0. 002,0 
52,63,69, 70 ,71,64,54,53,4,0. 002 , 0 . 002 , 0 
69,80, 86, 87,88,81 , 71 , 70, 4,0.002,0.002,0 

Neuman  Boundary  Condition  Cards 


None  specified,  so  no  cards  will  be  provided 


Dirichlet  Boundary  Condition  Cards 

1 . 110.5.2. 1 

36.55.5.2. 1 

The  first  node  of  the  first  Dirichlet  Boundary  is  node  1  and  the 
DBVAL  value  is  110  feet.  There  are  5  elements  on  this  face  that 
makes  this  Dirichlet  boundary  condition,  NDSAME  =  5.  NDINC  equals 
2  for  the  nodal  number  increase  for  the  consecutive  elements,  as 
in  3  -  1,  5-3,  7-5,  etc.  NMID  equals  1  for  the  node  numbering 

increase  between  successive  nodes. 

Free  Surface  Nodal  Cards 

11, 1,28, 18,45,35,62,52,79,69,96,36 

Here  the  total  mesh  is  allowed  to  move,  where  the  user  may  fix  the 
bottom  nodes  at  any  level  in  the  mesh. 

Mid-point  Free  Surface  Nodal  Cards 

17, 34, 51 , 68, 85 

Seepage  Face  Cards 

86,2,2 

The  seepage  is  prescribed  to  start  at  node  86,  NSFC  =  86,  and 
there  are  2  elements  for  this  seepage  face.  NSINC  is  equal  to  2 
as  there  is  a  two  nodal  increase  between  element  corner  nodes,  as 
is  88  -  86  and  90  -  88. 

Statement  and  End  Card 

END  OF  SILTY  SAND  DAM. 

3 , . 00 1 
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SILTY  SAND  DAM  . 

2,0.01 

12,5,0,6,50,2,1,0 
1,0,0, 1 
11,0,110,0 
18,200,0,1 
28,200,99,0 
35,400,0,1 
45,400,88,0 
52,6^0,0,1 
62,600,77,0 
69,800,0,1 
79,800,66,0 
86,1000,0,1 
96,1000,55,0 

1,12,18,19,20,13,3,2,4,0.002,0.002,0 

18,29,35,36,37,30,20,19,4,0.002,0.002,0 

35,46,52,53,54,47,37,36,4,0.002,0.002,0 

52,63,69,70,71,64,54,53,4,0.002,0.002,0 

69,80,86,87,88,81,71,70,4,0.002,0.002,0 

1.110.5.2.1 

86.55.5.2.1 

11,1,28,18,45,35,62,52,79,69,96,86 

17,34,51,68,85 

86,2,2 

END  OF  SILTY  SAND  DAM. 

3,0.01 
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.  aooo 

NOOI 

3  > 

1 1 

400 . OOOO 

Y  • 

31 

aooo 

NOOI 

4  1 

X  ■ 

400 . OOOO 

V  1 

13 

aooo 

NOOI 

43 

X  i 

400 . OOOO 

Y  * 

70 

4000 

NOOI 

41 

X  ■ 

400 . OOOO 

Y  • 

aa 

oooo 

NOOI 

■  3 

X  * 

aoo  oooo 

Y  ■ 

0 

oooo 

NOOI 

14 

X  * 

too  oooo 

V  i 

it 

4000 

NOOI 

•  • 

X  ■ 

aoo  oooo 

V  • 

30 

.  aooo 

NOOI 

•  • 

1  • 

aoo . oooo 

V  • 

4  a 

aooo 
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•  O 

X  • 

too  oooo 

V  • 

a  i 

aooo 

NODI 

•  3 

X  * 

too  oooo 

Y  • 

77 

oooo 

NOOI 

•  • 

X  i 

aoo  oooo 

V  ■ 

0 

oooo 

NOOI 

7  l 

1 1 

aoo . oooo 

Y  • 

1  3 

3  90 

NOOI 

73 

X  ■ 

aoo  oooo 

Y  • 

a  a 

4000 

NODI 

7  > 

XI 

aoo . oooo 

Y  • 

sa 

aooo 

NOOI 

77 

XI 

aoo . oooo 

V  ■ 

a  3 

aooo 

NOOI 

7  • 

1  • 

aoo  oooo 

Y  • 

aa 

oooo 

NOOC 

at 

X  ■ 

ooo . oooo 

Y  • 

0 

oooo 

NOOI 

•a 

ooo . oooo 

V  * 

1 1 

oooo 

NODI 

•  0 

ooo . oooo 

Y  • 

33 

oooo 

NOOC 

•  3 

ooo . oooo 

Y  • 

33 

oooo 

NOOI 

•  4 

Xi 

ooo  oooo 

Y  I 

44 

oooo 

NOOI 

•  a 

X  ■ 

ooo . oooo 

V  • 

aa 

oooo 

NEW  COORD  I NATES  OF  THE  FREE  SURFACE  LINE 


£ 


I 


V 

V 


i 

> 

I 


V. 


v*\ 


I 


NOOE  ■ 

71 

I 

too . OOOO 

V  ■ 

5371 

NODE  • 

1 2 

X 

•00  OOOO 

V  • 

•  1  .  77 • • 

NOOE  • 

AS 

X 

400  OOOO 

V  « 

S  2  3734 

NODE  • 

21 

X 

200  OOOO 

V  1 

10t  7  1  9  t 

NOOE  • 

1  1 

X 

0 . OOOO 

V  • 

1  I  0  OOOO 

CONVERGENCE  CHECK 

ERROR  ■ 

. 1 7 1 OE ■ 02 

NODE 

X  COORDINATE 

V  COORDINATE 

V  - FOTENT I AL 

X-VELOCITY 

Y- VELOCITY 

TOTAL  VeLO 

1 

0  OOOO 

0 . OOOO 

110  OOOO 

•  4314E04 

-.11  2SE -01 

. 43 1 4E -  04 

2 

0  OOOO 

1  1  . OOOO 

110  OOOO 

- . 43ISE  *04 

-  1  547E  -  1  5 

43 1 5E -  04 

3 

0  .  OOOO 

22 . OOOO 

no  oooo 

■  4304E -04 

. 425 1 E -07 

4304E -04 

4 

0  OOOO 

33  OOOO 

110. oooo 

•  4216E  04 

0 

4294E- 04 

s 

0  OOOO 

44  OOOO 

110  oooo 

-  .  4  2  7  9  E -04 

.  17S9E -0* 

4279E -04 

t 

0  OOOO 

SS  OOOO 

110  oooo 

-  4259E  04 

0 

. 82SSE -04 

7 

O  OOOO 

44  OOOO 

1  1  o  oooo 

-  4  2  2  7  E -04 

.  2920E -04 

. 4  2  2  7  E -04 

4 

0  OOOO 

77  OOOO 

1  io . oooo 

-  4 14 1 E  •  04 

0 

. 4 1 4 1 E -04 

1 

0  OOOO 

44  OOOO 

1 IO  oooo 

-  .  4  1  44  e • 04 

4499E -Of 

.414*9-04 

1  0 

0  OOOO 

41  OOOO 

no  oooo 

-  .  4  09  9  E -04 

losse - i b 

. 4099E -04 

1  1 

O  OOOO 

1  1  O  OOOO 

110  oooo 

-  .  4  0  4  5  E -  04 

7  5  7  7  E -04 

. 404SE -04 

1  2 

1 00  OOOO 

0  OOOO 

105  7701 

-  4  9  0  0 E -  04 

-  2  1  39E  -07 

. 4400E -04 

l  3 

1  00  OOOO 

21.1715 

105  7743 

-  459SE -04 

4574  E -  04 

. 459BE -04 

1  4 

1 OO  OOOO 

42  3430 

1047444 

-  . 45  42  E -04 

. 1 1 3 1 E -05 

4543E- 04 

t  5 

1 OO  OOOO 

434144 

104.4024 

- . 4540E • 04 

.  1 731 E -  05 

.45429-04 

1  1 

1  OO  OOOO 

44  4440 

1054243 

-  4530E -  04 

.  2  40  3  E -05 

. 4533E -04 

t  7 

1 OO  OOOO 

1  OS  4  S  7  5 

1 04  4434 

•  .  441  1 E -  04 

. 3223E -OS 

. 44979-04 

*  4 

200  OOOO 

0  OOOO 

lOI . 4005 

-  .  4  4  9  5  E -04 

. 1 1 4SE -07 

.  44BSE  -  04 

1  1 

200  OOOO 

i  o  ml 

1014014 

-  .  449  5E -  04 

.  3  4  7  3  E -04 

. 4  4  9  4  E -04 

20 

200  OOOO 

20 . 3430 

tOt  4044 

* . 449SE -04 

4 1 47E -04 

.  4  9  9  5  E -04 

2  t 

200  OOOO 

304144 

101  4092 

-  4  4 1 5  E -04 

.  1  1 S7E -05 

.  4494E -04 

22 

200  OOOO 

404440 

101 .4142 

•  441  SB  -  04 

1 4  3  OE *05 

.  44  9  4  E -04 

23 

200  OOOO 

SO  .  4575 

101 . 4244 

-  441 9E -04 

. 1909E -05 

.  4494E -  04 

24 

200  OOOO 

41  0241 

101  4354 

-  449  SE • 04 

243  1  E -  05 

.  849  BE  -  04 

24 

200  OOOO 

7  1  .  2004 

101.4440 

•  4  4  9  4  E -04 

2422E -05 

.  4902E -04 

24 

200  OOOO 

4  1  .372  1 

101.4422 

•  449  4  E -  04 

.  3  2  0  8  E  •  OS 

.  4  9  0  5  E -04 

27 

200  OOOO 

•  1  .4435 

101 .4740 

-  4404E -  04 

. 32SSE -05 

. 49 1 OE -04 

24 

200  OOOO 

10  1 .7151 

lOI  .  4  9  S  4 

-  .  4104E -  04 

. 3144E -05 

. 49 1 3E -04 

21 

300  OOOO 

O  OOOO 

94 . 4349 

-  .  9344E-04 

. 8744E -04 

. 9  3  5  4  E -04 

30 

300  OOOO 

14  4044 

944392 

-  . 9  3  5  4  E -04 

.  4  7  3  5  E -08 

.  93S5E -04 

3  l 

300  OOOO 

34.4177 

944414 

•  .  4  3  5  0  E • 04 

1 745E -05 

.43529-04 

32 

300  OOOO 

442244 

95 . 4730 

•  .  9343E -04 

.241 1 E -05 

.  934  BE  -  04 

33 

300  OOOO 

77 . 4344 

94 . 4024 

• . 9332E -04 

.  3  4  8  4  E -05 

.  9339E -04 

34 

300 . OOOO 

47.0442 

14 . 9404 

*.93149-04 

. 433 1 E -05 

.  432BE -04 

34 

400 . OOOO 

O  OOOO 

92 . 0447 

•  .  940 1 E -04 

-  4037E -04 

. 940 1 E -04 

31 

400 . OOOO 

4 . 2373 

92 . 0449 

- . 440 1 E -04 

.  492  1  E -04 

.94019-04 

37 

400  OOOO 

14 .4747 

92.0492 

-  .  4  4  00  E -  04 

.  9744E -04 

.  9400E -04 

34 

400 . OOOO 

27.7120 

92 . 0444 

-  .  9  7  4  4  E -04 

.  I47SE -04 

.  *7999-04 

31 

400  OOOO 

34.4414 

42.0424 

- . 979SE -04 

I957E-05 

. 9  7  4  7  E -04 

40 

400  OOOO 

44  1447 

92  073 1 

•  .  9742E  *04 

2  4  4  4  E *04 

.  979IE -04 

4  1 

400  OOOO 

45 . 4240 

92 . 0444 

•  .  4744E  *04 

.  2433E-09 

.  9792E -04 

42 

400  OOOO 

444414 

92.1 002 

•  .  9  7  4  3  E -04 

. 3423E -09 

.  9749E  *04 

43 

400 . OOOO 

73 . 4447 

92.1172 

-  .  97749  *04 

. 340BE -04 

.  9744E-04 

44 

400 . OOOO 

43  1341 

92  1342 

-  .  4772E  *04 

. 4  3  4  9  E -OS 

.  4742E-04 

44 

400  OOOO 

42 . 3734 

92  1474 

•  .  97SSE -04 

4474E -05 

.  97779-04 

44 

400  OOOO 

O . OOOO 

41  .  9474 

-  1043E-03 

.  2  1  44E -04 

1  0  4  3  E -03 

47 

400  OOOO 

17.4140 

44 . 9922 

-  1043E -03 

.  1  1  1 3E -04 

.  10*31-03 

4  1 

400  OOOO 

34 . 4300 

47  0044 

•  1  043E -03 

.  2223E-OS 

. 1043E -03 

41 

400  OOOO 

S  2 . 244  1 

47 . 0304 

•  1  042E -03 

. 33249-05 

.  10429 -03 

40 

BOO  OOOO 

44.4401 

47 . 0442 

-  1  040E -03 

4420E -05 

.  1 04 i e -os 

4  t 

500  OOOO 

47 . 075  1 

47  .  1073 

•  1  039E -03 

.  5S04E  -04 

. 1 040E -03 

52 

•OO  OOOO 

O  OOOO 

41.9114 

-  .  1  1 0  4  E -03 

-  5594E  *04 

1 1 OBE -03 

S3 

400  OOOO 

4  .  1777 

41.4129 

-  .  1  104E -03 

. 9374E -04 

.  1  1 OBE -03 

54 

800 

oooo 

1  4  3494 

81.9184 

- . 1 109E -03 

. 1 270E -04 

.  1  1 O  BE  -  03 

55 

BOO 

oooo 

24  5330 

4 1 . 4234 

•  .  1  104E -03 

. 19 10E -05 

. 1 104E -03 

58 

400 

oooo 

32.7107 

4 1  4325 

-.11 OBE -03 

2S34E -04 

.  1  1  O  BE  -  03 

57 

900 

oooo 

40  4444 

4  1  .  *44  1 

•  .  1  1O1E-03 

3 1 7 1 E  *05 

1  1 OBE -03 

54 

800 

oooo 

490481 

*  1  .  8444 

- . 1 104E -03 

. 3402E -05 

.  1  1058-03 

59 

400 

oooo 

57  2434 

4  1  .  *75  1 

-  .  1  1  03E -  03 

.  440BE -04 

.  1  1 04E -03 

40 

400 

oooo 

944214 

41.4944 

- . 1 1 02E -03 

.  40B3E  *05 

.  1  1031 -03 

9  1 

400 

oooo 

73  499  1 

41.7140 

•  .  1  10 1 E -  03 

5441E05 

.  1  1  0  3  E -03 

42 

900 

oooo 

8  1  .  77*4 

4  1 .740  1 

-  .  1  1  OOE -03 

. 429SE -04 

.  1  1  O  2  E -03 

43 

700 

oooo 

O  OOOO 

75.8434 

-  .  1  202E  *03 

.  32 1 7E -  09 

1 20  2  E • 0  3 

84 

700 

oooo 

14.  1314 

75 . 8444 

-  1 20  2  E -03 

.  1  49 1 E -04 

.  1  202E -03 

*5 

700 

oooo 

30  2824 

75  4455 

-  .  1 20 1 E  *  03 

. 29 1  BE -OS 

.  1  20  1  E -03 

4  4 

700 

oooo 

45 . 3442 

74 . 9932 

•  .  1 2  00  E -03 

4394E -05 

.  1  20 1 E -03 

87 

700 

oooo 

40  5244 

759322 

•  .  1  1  HE -03 

5  7  4  2  E • 05 

. 1 200E -03 

44 

700 

oooo 

75  4470 

75*427 

-  .  1  1  97E  *  03 

7  2  3  3  E -OS 

.  1  1  B9E -03 

89 

400 

oooo 

0  OOOO 

81.4922 

-  .  1 30 1 E -  03 

1 4 1 5E -07 

1301* -03 

70 

400 

oooo 

8  9937 

89  4*35 

-  .  1 30 1 E  *  03 

.  4  2  7  7  E *08 

. 1 30 1 E -03 

7  1 

400 

oooo 

13  4074 

41  9940 

-  .  1 3  OOE -03 

.  1  441E -05 

.  1  300E -03 

72 

400 

oooo 

20  44  11 

49  9049 

-  1 299E -03 

24  4  1  E  •  OS 

.  1  294E -03 

73 

800 

oooo 

27  4149 

498141 

-.12991*03 

. 387 1 E -05 

1  2  *  4  E -03 

74 

400 

oooo 

34  7444 

41  4274 

-  .  1 21  BE -03 

4  0  2  4  E -  05 

.  1  297E -03 

75 

400 

oooo 

4  1  7  2  2  3, 

498431 

•  .  1294E -03 

. 543 1 E -05 

. 1 295E -03 

74 

400 

oooo 

44  97  40' 

49  4404 

-  .  1  242E  -03 

S4SSE -05 

.  1  293E -03 

77 

400 

oooo 

SB  4297 

99.9410 

•  .  1  28BE -03 

.  7043E -  05 

1  29  OE -03 

7* 

400 

oooo 

82.4434 

49  7030 

-  1 244E -03 

4570E  -05 

1  247E -  03 

79 

400 

oooo 

445371 

41  7274 

-  1  242* -  03 

8  870E  -05 

1  245E -03 

40 

900 

oooo 

0  OOOO 

42  9454 

- . 1 4S1E -03 

•11 OBE -09 

.  1  459E -  03 

8  1 

900 

oooo 

12  4537 

42  4425 

-  1 4  SOI  -  03 

.  1  723E • 05 

1  480E -03 

*2 

•  00 

oooo 

24.9074 

427121 

-  .  1 4  4  0  E -03 

.  3  4  7  BE  -  05 

1 44 1 E -03 

43 

100 

oooo 

37  34  11 

42  7474 

-  1442E  -03 

5  2  9  O E -  05 

1  443E -03 

44 

900 

oooo 

494144 

42  7472 

•  1  444E  •  03 

7 1  *  BE  -  05 

.  1  4  4  4  E -03 

45 

400 

oooo 

42  2444 

42.4434 

-  1 4  44E -  03 

•  4  3  7  B -05 

f 444E-03 

44 

1  OOO 

oooo 

0  OOOO 

54  OOOO 

•  .  1  4  1  SE -  03 

-  1747E-08 

1  9  1  5E -03 

4  7 

1  OOO 

oooo 

5 . 4000 

SB  OOOO 

- . 1 4 1 8E • 03 

10S4E  -  IS 

1  9  1  BE  -  03 

94 

1  OOO 

oooo 

1 1  OOOO 

55  OOOO 

-  1 4 1 7E • 03 

502 1 E -08 

.  1  4  1  7E -  03 

49 

1  OOO 

oooo 

1 4  5000 

55  OOOO 

-  1  420E -  03 

-  1  OS4E -  15 

1  420E -03 

•  0 

1  OOO 

oooo 

22  OOOO 

54  OOOO 

-  1  423E -  03 

1  034E -05 

.  1  S23E -03 

9  1 

1  OOO 

oooo 

27  4000 

44  OOOO 

•  1  421E  •  03 

O 

1 429E -03 

92 

1  OOO 

oooo 

33  OOOO 

55  OOOO 

-  1  434E -  03 

.  1  4  2  7  E -05 

.  1  434E -03 

93 

1  OOO 

oooo 

34  5000 

55  OOOO 

-  1  442E  -03 

0 

.  1  442E -03 

94 

1  OOO 

oooo 

44  OOOO 

54 . OOOO 

-  1  4  4  4  E  -03 

23 1 9E -05 

.  1  4  5  OE -03 

95 

1  OOO 

oooo 

44  5000 

54 . OOOO 

-  .  1  81  1  E  -  03 

0 

1  44  1  E -  03 

94 

rooo 

oooo 

SI  OOOO 

SS . OOOO 

•  .  1  444E  -03 

.  3  3  2  4  E -OS 

1  9448  -  03 

S  E  E  F  A  C  E  FACE  CONFUTATION  RESULTS 


N  4 

49  NE  *  44 

AVE 

VEL  • 

1  4  1  40E -03 

OS  * 

.  1  1 OOOE  *02 

DO  * 

.  1  7778E -02 

N4 

TOTAL  0 I SCHARCE 
»  44  ME  *  90 

AYE 

1 7774E 
VEL  * 

•02 

1  4  2 0 2 E -  03 

DS  * 

1  1 OOO E  * 0 2 

DO  * 

.  1  7422E -02 

TOTAL  DISCHARGE 

* 

-02 

■NO  OF  SI 

OININ  •  3 

LTV  9. 

AND  OAM 

FINN  • 

0 1 000 

A- 19 


FRES  COEF 


T 
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APPENDIX  B— 1 


SAMPLE  INPUT 


AND  OUTPUT 


SP  -  25  FEET 


DEPTH  -  20  FEET  K  -  .2.  DEPTH  -58'. 


g 


I 


i; 

s.. 


R 


zv 


i 


$ 


l 


§ 


* 


i 


I 


RUN  ADS 2 

2,0.0001 

93,30,0,30,20,6,5,0 

I, -340,0,1 
5,-340,34,1 

II, -340,78.0,0 
18,-200,0,1 
22,-200,34,1 
28,-200,78.0,0 
35,-80,0,1 
39,-80,34,1 
45,-80,77.0,0 
52,0,0,1 
56,0,34,1 
62,0,74.0,0 
69,50,0,1 

73.50.34.1 
79,50,71.0,0 
86,58.75,0,1 

90.58.75.34.1 
96,59.0,58.0,0 
103,62.0,0,1 
107,62.0,34,1 
113,61.0,58.0,0 
120,66.0,0,1 
124,66.0,34,1 
130,65.0,60.0,0 
137,70.0,0,1 
141,70.0,34,1 
147,69.0,60.0,0 
154,74.0,0,1 
158,74.0,34,1 
164,75.0,60.0,0 
171,78.0,0,1 
175,78.0,34,1 
181,79.0,60.0,0 
188,82.0,0,1 
192,82.0,34,1 
198,83.0,58.0,0 
205,86.0,0,1 
209,86.0,34,1 
215,85.0,58.0,0 
222,90.0,0,1 
226,90.0,34,1 
232,89.0,60.0,0 
239,94.0,0,1 
243,94.0,34,1 
249,93.0,60.0,0 
256,98.0,0,1 
260,98.0,34,1 
266,99.0,60.0,0 
273,102.0,0,1 
277,102.0,34,1 
283,103.0,60.0,0 
290,106.0,0,1 
294,106.0,34,1 
300,107.0,58.0,0 
307,110.0,0,1 
311,110.0,34,1 
317,109.0,58.0,0 


324,114.0,0,1 

328,114.0,34,1 

334,113.0,60.0,0 

341,118.0,0,1 

345,118.0,34,1 

351,117.0,60.0,0 

358,122.0,0,1 

362,122.0,34,1 

368,123.0,60.0,0 

375,126.0,0,1 

379,126.0,34,1 

385,127.0,60.0,0 

392,130.0,0,1 

396,130.0,34,1 

402,131.0,58.0,0 

409,134.0,0,1 

413,134.0,34,1 

419,133.0,58.0,0 

426,138.0,0,1 

430,138.0,34,1 

436,137.0,60.0,0 

443,142.0,0,1 

447,142.0,34,1 

453,141.0,60.0,0 

460,146.0,0,1 

464,146.0,34,1 

470,147.0,60.0,0 

477,150.0,0,1 

481,150.0,34,1 

487,151.0,60.0,0 

494,154.0,0,1 

498,154.0,34,1 

504,155.0,58.0,0 

511,156.0,0,1 

515,156.0,34,1 

521,156.0,58.0,0 

1,12,18,19,20,13,3,2,4,0.2,0 

18,29,35,36,37,30,20,19,4,0. 

35,46,52,53,54,47,37,36,4,0. 

52,63,69,70,71,64,54,53,4,0. 

69,80,86,87,88,81,71,70,4,0. 

86,97,103,104,105,98,88,87,4 

103,114,120,121,122,115,105, 

120,131,137,138,139,132,122, 

137,148,154,155,156,149,139, 

154,165,171,172,173,166,156, 

171 , 182, 188,189, 190, 183, 173, 

188,199,205,206,207,200,190, 

205,216,222,223,224,217,207, 

222 , 233 ,239,240,241,234,224, 

239,250,256,257,258,251,241, 
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409 , 420 , 426 , 427 , 428 ,421,411,410,4,0.2,0.2,0.0 
426,437,443,444,445,438,428,427,4,0.2,0.2,0.0 
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